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EXECUTIVE  SUMMARY 


The  removal  of  nitric  oxide  (NO)  from  air  streams  using  pulsed  streamer  corona 
discharge,  a  type  of  non-thermal  plasma,  is  considered  in  the  present  study.  Pulsed  streamer 
corona  technology  is  based  upon  the  application  of  an  electric  discharge  to  produce  highly 
reactive  free  radicals  in  a  gaseous  medium.  In  this  study,  a  pulsed  streamer  corona  discharge 
reactor  was  operated  to  study  the  removal  of  NO-containing  gases.  The  reactor  was  designed, 
constructed,  and  preliminary  analysis  and  experiments  performed  during  previously  funded 
phases  of  this  project.  Most  of  the  details  concerning  these  aspects  were  reported  in  the  final 
report  documents  for  those  phases.  In  the  present  study  physical  and  chemical  characterization  of 
the  pulsed  corona  reactor  were  performed.  The  physical  characteristics  of  the  discharge,  such  as 
the  rise  time,  peak  voltage,  current,  power,  and  pulse  width  of  the  voltage  pulse  were  measured. 
NO  removal  from  feed  gases  of  various  compositions  including  pure  nitrogen,  and  dry  and  humid 
air,  with  and  without  added  ethylene  and  carbon  monoxide,  was  studied.  In  an  atmosphere  of 
nitrogen,  the  main  reaction  mechanism  for  the  removal  of  NO  was  found  to  be  through  reduction 
of  NO  to  molecular  nitrogen,  while  in  the  case  of  dry  air,  NO  was  oxidized  to  NO..  In  humid 
air,  NO.  removal  occurs  through  oxidation  pathways  to  nitric  acid.  The  addition  of 
hydrocarbons,  namely  ethylene,  significantly  affected  NO  removal. 

A  kinetic  reactor  model  was  also  developed  to  characterize  the  chemical  reactions  taking 
place  in  the  pulsed  streamer  corona  reactor.  The  kinetic  rate  constants  for  the  production  of 
highly  reactive  free  radicals  such  as  O,  N,  and  OH  were  obtained  by  fitting  the  model  to  the  NO 
removal  experimental  data.  The  model  was  then  used  to  predict  the  output  concentrations  of 
various  other  species,  including  O3,  N.O,  and  a  number  of  organic  compounds,  in  order  to  test 
the  model  and  to  guide  the  experimental  analysis  of  byproduct  characterization. 

The  major  findings  of  the  present  work  include: 

1)  Task  1  (Section  7.1.4  RFP)  -  Gas  Phase  Pulsed  Corona  Experiments. 

An  extensive  series  of  experiments  to  consider  the  effect  of  the  feed  gas  composition  on 
the  pulsed  corona  reactor  performance  for  a  range  of  reactor  residence  times  and  applied  pulsed 
voltages  was  conducted.  Experiments  with  NO  in  dry  nitrogen,  NO  in  dry  air,  NO  in  humid  air. 
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NO  and  ethylene  in  dry  and  humid  air,  and  NO,  ethylene,  and  CO  in  dry  and  humid  air  were 
performed  with  reactor  residence  times  varying  jfrom  approximately  5  seconds  to  60  seconds  and 
applied  pulsed  voltages  ranging  from  30  kV  to  50  kV.  A  minimum  of  two  repeat  runs  for  each 
set  of  conditions  were  performed.  Specific  results  from  these  runs  indicate: 

a)  In  pure  nitrogen  gas,  NO  was  chemically  reduced  to  Nj.  About  99%  removal  of  NO„ 
was  observed  at  50  kV  pulsed  voltage  for  a  gas  residence  time  of  44  seconds  in  pure  nitrogen. 

NO  removal  in  pure  nitrogen  occured  at  a  lower  energy  density  than  that  reported  in  the  literature 
for  dielectric  barrier  discharge  reactors  and  some  corona  discharge  reactors  (Penetrante  et  al., 
1996).  In  addition,  the  pure  nitrogen  experiments  tend  to  support  the  hypothesis  that  the  NO 
removal  is  dependent  not  only  on  the  power  density  to  the  reactor,  but  also  upon  the  absolute 
value  of  the  applied  electric  field.  Wolfetal.  (1996),  Mizuno  et  al.  (1995),  and  Masudaetal. 
(1990)  also  found  that  the  field  strength  affects  the  NO  removal  in  pulsed  corona  discharge  and 
dielectric  barrier  reactors  independent  of  the  power  density,  in  contradiction  to  the  theoretical 
predictions  and  claims  of  Penetrante  et  al.  (1996).  Our  data  for  NO  removal  in  nitrogen  matches 
quite  closely  with  that  reported  by  Wolf  et  al.  (1996)  for  their  work  at  high  electric  fields  with 
dielectric  barrier  discharge  reactors. 

b)  All  NO  removal  rates  were  monotonic  functions  of  the  reactor  residence  time,  and 
showed  increasing  amounts  of  NO  removal  with  higher  voltages.  NO  removal  from  dry  air  in 
the  presence  of  ethylene  as  a  function  of  reactor  residence  time  showed  the  fastest  removal  rate 
of  all  the  conditions  studied.  The  addition  of  water  vapor  to  the  feed  gas  (saturated  at  room 
temperature)  did  not  enhance  the  NO  removal.  The  concentration  profiles  of  NO  versus 
residence  time  for  the  various  applied  voltages  in  humid  air,  in  the  presence  of  500  ppm  CO  in 
humid  air,  and  in  the  presence  of  500  ppm  of  ethylene  in  humid  air  were  almost  the  same. 

c)  The  oxidation  of  NO  leads  to  the  formation  of  NO2,  and  therefore  many  of  the 
experiments  showed  an  increase  in  NO2  concentration  with  residence  time  followed  by  a 
decrease  at  very  long  residence  time,  as  expected  for  a  sequence  of  chemical  reactions  occurring 
in  series.  In  addition,  generally  more  NO2  was  formed  at  the  higher  voltages  in  the  presence  of 
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oxygen  but  not  in  pure  nitrogen.  The  addition  of  oxygen  in  the  form  of  dry  air  increased  the  rate 
of  NO  removal,  however,  it  also  increased  the  rate  of  formation  of  NOj  as  expected  by  the 
chemieal  kinetics.  In  dry  air,  NO  was  oxidized  to  form  NO2  and  N2O.  NO2  formation  was 
significantly  reduced  in  the  presence  of  water  vapor.  The  addition  of  ethylene  to  the  dry  air  feed 
gas  enhanced  NO2  removal.  In  humid  air,  however,  the  addition  of  500  ppm  CO  and  500  ppm 
ethylene  had  almost  no  effect  on  NO2  removal  compared  to  the  case  of  humid  air  alone. 

2)  Task  2  (Section  7.1.2  RFP)  and  Task  4  (Section  7.1.6  RFP)  -  Reactor  and  Kinetic  Modeling. 

Chemical  reaction  modeling  was  performed  using  the  CHEMKIN  software  package  with 
details  described  in  the  Modeling  section  of  this  report.  The  specific  results  from  the  model  and 
model/data  comparison  include: 

a)  In  general,  the  model  profiles  of  NO  and  NO2  matched  the  experimental  data  fairly 
well,  and  the  model  indicated  a  small  amount  of  N2O  formation  which  was  also  observed 
experimentally.  The  overall  reaction  rate  constant  for  nitrogen  dissociation  was  determined  by 
fitting  the  model  with  the  NO  experimental  data.  The  experiments,  as  well  as  the  modeling 
results,  follow  the  general  trends  indicated  in  the  literature  on  non-thermal  plasma  treatment 
methods. 

b)  The  rate  constants  for  nitrogen  and  oxygen  dissociation  were  determined  from  the  NO 
removal  and  the  NO2  formation  data.  The  reaction  rate  constant  for  nitrogen  dissociation  was 
different  from  the  value  obtained  in  pure  nitrogen;  this  is  also  expected  based  upon  observations 
reported  in  the  literature  that  electronegative  gases  affect  the  electron  energy  distribution  and 
thus  the  chemical  reaction  rate  constants.  The  pseudo  first-order  rate  constant  for  nitrogen 
dissociation  was  0.0475  s''  in  an  atmosphere  of  nitrogen  at  50  kV,  but  in  dry  air  the  rate  constant 
was  lowered  to  0.03  s  '.  Hence,  as  expected  from  the  literature,  the  rate  constant  depends  upon 
the  composition  of  the  gas.  The  differences  in  the  rate  constants  for  nitrogen  and  oxygen 
dissociation  in  dry  air  were  functions  of  the  applied  voltage.  At  low  voltage  (30  kV)  the  rate 
constants  were  almost  equal  (approximately  .01  s'')  and  at  the  high  voltage  (50  kV)  the  rate 
constants  were  0.025  s''  for  nitrogen  dissociation  and  0.25  s''  for  oxygen  dissociation.  This 


vi 


factor  often  difference  is  significantly  smaller  than  theoretical  studies  (Penetrante  et  al.,  1995) 
where  it  was  predicted  that  the  G-values  for  oxygen  and  nitrogen  dissociation  would  differ  by 
about  600  times.  It  should  also  be  noted  that  these  theoretical  studies  also  do  not  predict  a  field 
dependence  as  seen  in  our  work. 

c)  The  measured  N2O  matched  well  with  the  predictions  of  the  model.  This  result  is 
significant  since  this  is  an  independent  check  of  the  model.  The  amount  of  ozone  produced  from 
dry  air,  as  measured  using  a  KJ  titration  method,  was  about  a  factor  of  three  higher  than  that 
predicted  by  the  model.  Further  refinement  of  the  model  and  additional  experiments  to  measure 
ozone  concentration  is  continuing. 

3)  Task  3  (Section  7. 1 .5  RFP)  -  Byproduct  Identification. 

a)  Complete  removal  of  ethylene  was  observed  for  the  case  of  humid  air  with  NO.  The 
exhaust  gas  from  the  reactor  was  sparged  through  a  liquid  solution  and  this  liquid  was  analyzed 
by  HPLC.  Byproduct  analysis  by  HPLC  indicated  that  less  than  1%  of  the  carbon  from  500  ppm 
ethylene  leads  to  the  formation  of  formic  acid.  Two  additional  peaks  of  similar  magnitude  to  the 
formic  acid  peak  were  observed  on  the  HPLC.  Through  a  process  of  elimination,  it  was  shown 
that  these  peaks  were  not  due  to  acetic  acid,  formaldehyde,  ethylene  glycol,  glycolaldehyde, 
ethanol,  or  methanol.  Exact  identification  of  these  peaks  was  not  successful;  however,  based  on 
the  magnitude  of  the  formic  acid  peak,  it  is  expected  that  these  peaks  represent  at  most  a  few 
percent  of  the  initial  total  carbon  from  the  ethylene  degradation.  Work  is  continuing  to  make 
further  identificaiton  of  these  peaks. 

b)  The  liquid  samples  were  also  analyzed  by  ion  chromatography  where  nitrate  and  nitrite 
could  be  quantified.  This  analytical  method,  which  is  very  accurate  and  has  a  detection  limit  in 
the  ppb  range,  showed  no  nitrite  and  a  small  amount  of  nitrate.  The  amount  of  nitrate  measured 
by  this  technique  corresponds  to  about  5%  of  the  nitrogen  from  the  feed  NO  converted  to  nitric 
acid.  It  thus  appears  that  a  balance  on  nitrogen  was  not  obtained.  Of  about  100  ppm  NO  in  the 
initial  feed  gas  of  a  humidified  gas  stream,  only  a  few  percent  of  the  nitrogen  was  detected  as 
nitric  acid.  It  is  possible  that  some  of  the  acid  formed  could  have  been  deposited  on  the  walls  of 
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the  reactor,  however  it  was  difficult  to  assess  this  quantitatively.  Recommendations  for  future 
work  include  the  analysis  of  aerosol  formation. 

c)  In  dry  air,  over  95%  of  the  500  ppm  ethylene  fed  into  the  reactor  was  decomposed  at 
50  kV  pulsed  voltage  and  44  seconds  residence  time.  Approximately  130  ppm  of  CO2  was 
produced  during  these  conditions  as  measured  by  the  GC/MS.  This  would  imply  that 
approximately  13%  of  the  carbon  fed  to  the  reactor  left  the  reactor  as  carbon  dioxide. 
Measurements  by  HPLC  found  no  formic  acid  under  these  conditions.  Further  work  is  needed  to 
determine  the  fate  of  the  remaining  carbon. 

4)  Task  5  (Section  7.1.3  RFP)  -  Liquid  Reactor  System  Specification. 

A  liquid  phase  reactor  has  been  designed  and  constructed.  This  reactor  consists  of  a  1 
liter  jacketed  glass  vessel  equipped  with  ports  to  allow  the  pulsed  voltage  electrode  assembly  to 
be  inserted.  Modification  of  the  high  voltage  electrode  will  allow  rapid  changes  in  the  discharge 
tips.  The  top  of  the  reactor  will  allow  for  the  introduction  and  removal  of  samples  as  well  as  for 
electrical  connections  to  the  ground  electrode.  The  basic  features  of  this  reactor  are  similar  to 
those  used  in  previous  work  by  our  group  (Sharma  et  al.,  1993;  Joshi  et  al.,  1995),  but  the  reactor 
is  somewhat  larger  and  more  easily  operated  and  modified  than  previous  designs. 
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SECTION  I 
INTRODUCTION 


Nitrogen  oxides  released  into  the  atmosphere  arise  primarily  from  the  combustion  of  fossil 
fuels  from  both  stationary  and  mobile  sources.  In  1992,  U.S.  emissions  of  anthropogenic  nitrogen 
oxides  were  estimated  to  be  23  million  tons.  Of  this,  45%  was  from  mobile  sources  and  the 
remainder  came  from  stationary  sources  (Ozkan  et  ah,  1995).  The  seven  oxides  of  nitrogen  that  may 
be  present  in  ambient  air  are  NO,  NOj,  Np,  NO3,  NPj,  N2O4  and  NPj .  Of  these  nitrogen  oxides, 
the  only  ones  subject  to  regulatory  control  are  NO  and  NO2,  and  these  are  together  referred  to  as 
NO,.  Although  N2O  is  not  considered  to  be  an  air  pollutant,  it  contributes  to  stratospheric  ozone 
destruction  and  it  is  a  greenhouse  gas  (U.S.  EPA,  1993). 

Nitric  oxide  (NO)  is  an  odorless  gas  that  is  only  slightly  soluble  in  water  (0.006  g/100  g  of 
water  at  24  °C  and  1  atm  pressure).  It  is  toxic  and  makes  up  90-95%  of  the  NO,  emissions  from 
fossil  fuel  combustion.  Nitrogen  dioxide  (NO^)  is  a  reddish  brown  gas  with  a  characteristic  pungent 
odor.  It  is  corrosive  and  highly  oxidizing.  It  is  formed  rapidly  in  the  atmosphere  from  the  oxidation 
of  NO,  and  it  is  a  toxic  agent  which  can  negatively  affect  human  health.  Nitrous  oxide  (N2O),  or 
laughing  gas,  is  a  colorless  gas  with  a  slight  odor  at  high  concentrations  (U.S.  EPA,  1993). 

NO  and  NOj  are  highly  toxic  if  inhaled  in  high  concentrations.  NO  entering  the  body  forms 
nitrite  which  oxidizes  the  iron  in  hemoglobin,  rendering  it  ineffective  as  an  oxygen  carrier.  In  1986, 
the  World  Health  Organization  recommended  that  the  exposure  of  humans  to  NO,  should  not  exceed 
200  ppbv  for  more  than  one  hour  a  month.  Long-term  exposure  at  50  ppbv  can  cause  pulmonary 
damage  and  increased  airway  resistance  in  healthy  individuals  (U.S  EPA,  1993). 

NO,  gases  are  one  of  the  main  contributors  to  the  formation  of  photochemical  smog.  This 
smog  is  primarily  due  to  the  reaction  of  NO,  and  hydrocarbons  which,  at  high  concentrations, 
contribute  to  a  reduction  in  visibility.  Under  the  influence  of  UV  light  these  compounds  produce  a 
variety  of  harmful  compounds  such  as  peroxyacetyl  nitrate  (PAN)  and  ozone.  Ozone  is  a  highly 
toxic  substance  which  leads  to  breathing  problems  above  concentrations  of  500  ppbv.  NO, 
emissions  are  one  of  the  primary  contributors  to  acid  rain,  which  is  associated  with  a  number  of 
effects  including  acidification  of  lakes  and  streams,  and  accelerated  corrosion  of  buildings  and 
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monuments  (Leslie  et  al.,  1992). 

The  two  main  mechanisms  that  are  primarily  responsible  for  NO^  formation  during 
combustion  process  are  thermal  NO^  and  fuel  NO^.  Thermal  NO^  results  from  the  oxidation  of  N2 
contained  in  combustion  air,  whereas  fuel  NO,  originates  from  the  oxidation  of  chemically  bound 
nitrogen  in  the  fuel.  Abatement  of  NO,  emissions  can  be  achieved  by  (1)  modifying  the  combustion 
process,  such  as  using  low-NO,  burners,  gas  recirculation,  and  staged  combustion,  and  (2)  using 

post-combustion  removal  methods  (Chang  et  al.,  1992). 

Among  the  post-combustion  control  methods,  selective  catalytic  reduction  (SCR)  is  the  most 
advanced  technology  currently  used  to  remove  NO,  from  gas  streams.  One  disadvantage  with  the 
SCR  process  is  that  the  particulate  matter  commonly  found  in  flue  gas  streams  generated  by  the 
combustion  of  fossil  fuels  tends  to  "poison"  catalysts,  and  therefore  lower  SCR  performance. 
Another  post-combustion  removal  method  for  NO,  is  selective  non-catalytic  reduction  (SNCR). 
SNCR  relies  on  the  injection  of  ammonia  into  the  gas  stream  to  chemically  reduce  NO  to  N2  and  H2O 
in  the  presence  of  O2.  However,  SNCR  is  effective  only  between  900  and  1 100  °C  (Kokkinos  et  al., 

1991). 

One  alternative  to  SCR  and  SNCR  processes  that  has  been  studied  since  the  1 970's  is  the  gas 
phase  plasma  oxidation  of  NO,  using  an  electric  discharge.  Electric  discharges  initiate  electron-gas 
reactions  that  in  turn  produce  highly  reactive  radicals.  The  application  of  an  electric  discharge  to  a 
gas  can  lead  to  the  creation  of  a  thermal  or  a  non-thermal  plasma. 

A  thermal  plasma  is  a  "hot"  or  "equilibrium"  plasma,  which  is  characterized  by  a  high  gas 
temperature  and  approximately  equal  gas  and  electron  temperatures.  Typical  examples  are  those 
produced  in  arcs  and  plasma  torches.  A  "non-thermal  plasma",  as  the  name  implies,  is  a  plasma  in 
which  the  electron  temperature  is  considerably  higher  than  that  of  the  components  of  the  ambient 
gas.  The  electrons  produced  are  short-lived  and  rarely  collide  with  a  pollutant  molecule,  but  they 
undergo  many  collisions  with  the  predominant  gas  molecules,  producing  radicals  that  in  turn  lead 
to  the  removal  of  the  toxic  molecules  (Penetrante,  1993).  Hence,  the  essence  of  the  non-thermal 
plasma  technique  is  the  efficient  use  of  electrical  energy  to  generate  gas  phase  radicals  by  electron 

impact  excitation,  ionization,  and  dissociation. 

The  different  types  of  non-thermal  plasmas  that  have  been  investigated  are; 
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1.  Electron  Beam 

2.  Dielectric  Barrier  Discharge 

3.  Corona  Discharge. 

The  electron  beam  process  involves  the  direct  irradiation  of  gases  by  an  external  source  beam 
of  highly  energetic  electrons  to  produce  active  radicals  and  atoms,  which  react  with  NO^  to  form 
nitric  acid  as  well  as  other  products.  The  use  of  electron  beams  to  prevent  pollution  was  initiated 
in  the  1970's  in  Japan  by  the  Japan  Atomic  Energy  Research  Institute  and  the  Ebara  Corporation,  and 
pilot  scale  studies  are  currently  being  conducted  in  several  countries  such  as  the  U.S.,  Germany,  and 
Poland  (Frank  et  al.,  1993). 

In  a  dielectric  barrier  discharge,  the  plasma  is  sustained  between  one  or  more  dielectric 
surfaces  backed  by  a  conductor.  A  repetitively  pulsed  or  AC  voltage  is  applied  to  the  electrodes, 
resulting  in  a  filamentary  periodic  pulsed  plasma  (Chang  et  al.,  1992).  The  dielectric  barrier 
discharge  is  currently  used  in  industry  for  ozone  synthesis  (Dhali  et  al.,  1991;  Kogelschatz  and 
Eliasson,  1995). 

Corona  discharges  such  as  pulsed  streamer  corona,  AC  corona,  and  DC  corona,  have  been 
studied  for  the  treatment  of  flue  gas,  VOC's,  and  for  the  synthesis  of  gases  such  as  ozone.  In  the 
present  study,  the  removal  of  NO  from  dry  and  humid  air  streams  is  investigated  using  positive 
pulsed  streamer  corona.  The  effect  on  NO  removal  of  additives  such  as  carbon  monoxide  and 
ethylene  is  also  considered.  The  exhaust  gases  coming  out  of  jet  engines  (FI  01  engine  with  JP4 
fuel)  have  four  compounds  (ethylene,  acetylene,  propene,  and  formaldehyde)  accounting  for  27% 
of  volatile  organic  emissions  at  idle  power  (Spicer  et  al.,  1992).  Ethylene  is  used  in  the  present  study 
as  a  model  hydrocarbon  to  consider  the  effect  of  hydrocarbons  on  NO^  chemistry. 

A.  RESEARCH  OBJECTIVES 

The  specific  goals  of  this  research,  as  outlined  in  the  statement  of  work  for  this  contract,  include: 
Gas-Phase  Pulsed  Corona  Experiments. 

1)  Task  1.  Investigation  of  the  intermediate  and  final  by-products  formed  in  the  destruction  of 
combustion  source  exhausts.  Use  of  digital  storage  oscilloscope  or  high  speed  computer  data 


acquisition  system  to  measure  gas  phase  pulsed  corona  discharge  characteristics. 

Gas-Phase  Chemical  Kinetic  Analysis  and  By-Product  Correlation. 

2)  Task  2.  Develop  or  modify,  as  needed,  plasma  chemistry  codes  for  use  in  modeling  chemical 
kinetic  reactions  and  for  comparison  of  experimental  results  with  code  predictions. 

3)  Task  3.  Identify  and  estimate  quantities  formed  of  by-products  generated  during  the  experiments 

of  Task  1. 

4)  Task  4.  Application  of  codes  developed  during  the  performance  of  Task  2  to  simulate 
representative  experimental  data  generated  during  the  performance  of  Tasks  1  and  2. 

Liquid-Phase  Pulsed  Corona  Experiments. 

5)  Task  5.  Develop  and  design  a  pulsed  corona  discharge  reactor  for  use  in  aqueous  solutions. 

B.  BACKGROUND 

This  section  reviews  previous  research  on  pollutant  removal  from  waste  streams  using  non- 
thermal  plasma  techniques  such  as  electron  beams,  dielectric  barrier  discharge,  and  corona  discharge. 
Additional  details  are  reported  in  the  final  report  for  Phase  I  of  this  project.  Details  on  the  operation 
and  conduction  of  liquid  phase  experiments  will  be  reported  in  the  next  report  due  in  May  1997  as 
required  by  the  contract. 

Electron  Beams 

Electron  beam  treatment  to  prevent  pollution  was  first  initiated  in  Japan  in  1970  by  the  Japan 
Atomic  Energy  Research  Institute  in  collaboration  with  Ebara  Corporation.  Several  pilot  plant 
studies  of  electron  beam  process  applications  to  coal-fired  powered  plants  have  been  made  in  Europe 
and  in  the  U.S.  This  process  is  particularly  attractive  for  the  treatment  of  large  volumes  of  gaseous 
effluents  such  as  those  coming  out  of  utility  power  plants  (Frank  et  al,  1993). 

In  the  electron  beam  process,  irradiation  of  the  flue  gas  with  energetic  (300-800  keV)  electrons 
initiates  electron-gas  reactions  to  produce  radicals  that  have  high  enough  energ>'  to  oxidize  NO.,  and 
SO,  to  nitric  acid  and  sulfuric  acid,  respectively.  In  the  presence  of  ammonia  (NH3),  these  acids  are 
converted  into  ammonium  sulfate  and  ammonium  nitrate.  These  salts  are  later  recovered  as  a  dry 
powder  using  a  conventional  particle  collector  (Frank  et  al.,  1993,  Matzing,  1991). 
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The  reaction  mechanism  involved  in  the  removal  of  NO^  and  SO2  in  electron  beam  treatment  was 
investigated  by  Tokunaga  et  al.  (1984).  Reactions  of  NO  and  SO2  were  studied  under  the  irradiation 
of  electron  beams  using  various  gas  compositions  ranging  from  a  simple  system  (nitrogen-NO)  to 
a  more  complex  system  (nitrogen-oxygen-water  vapor-NO-S02)  which  is  close  to  the  composition 
of  an  actual  flue  gas.  For  the  NO-nitrogen  gas  mixture  at  a  gas  temperature  of  100  “C,  a  dose  of  3 
Mrad,  and  an  initial  NO  concentration  of  250  ppm,  10%  of  the  NO  was  converted  into  NO2  and  3% 
of  the  NO  was  converted  to  N2O.  The  remaining  NO  was  anticipated  to  be  decomposed  mainly  to 
nitrogen  and  oxygen.  In  the  case  of  a  N02-nitrogen  mixture  at  a  dose  of  6  Mrad,  60%  of  the  NO2 
formed  NO  and  N2O,  and  40%  of  the  NO2  was  decomposed  to  nitrogen  and  oxygen.  For  different 
NO  and  NO2  concentrations  in  a  gas  containing  3%  oxygen  and  97%  nitrogen  at  100  "C,  the  final 
concentrations  of  NO  and  NO2  were  represented  by  the  following  empirical  equation: 

M^=Qmippm-^  (1) 

[NO^f 

The  total  NO  concentration  was  observed  to  increase  with  dose  rate  for  initially  low 
concentrations  of  NO  (50-100  ppm),  and  was  not  observed  to  follow  the  above  empirical  relation. 
The  presence  of  water  vapor  (2.2%  to  12.9%)  enhanced  the  NO  removal  and  the  formation  of  HNO3 
at  a  gas  temperature  of  120  °C.  The  addition  of  additives  such  as  carbon  monoxide  (0-1.13  %)  and 
ammonia  (0-2000  ppm)  were  observed  to  have  a  positive  effect  on  the  removal  of  NO.  For  a  1 . 1 3% 
addition  of  carbon  monoxide,  NO  removal  was  enhanced  by  50%  at  1  Mrad,  and  an  increase  of  14% 
was  found  in  the  NO  removal  efficiency  with  the  addition  of  2000  ppm  of  NH3. 

A  kinetic  model  for  the  radiation-induced  oxidation  of  nitrogen  oxides  to  include  temperature 
dependence  and  the  presence  of  sulfur  dioxide  was  developed  by  Busi  et  al.  (1987).  An  increase  in 
temperature  from  ambient  to  120  "C  decreased  the  removal  efficiency  of  NO.  Matzing  (1991) 
considered  detailed  theoretical  and  model  results  for  the  electron  beam  dry  scrubbing  process 
(EBDS).  The  major  radical  source  was  shown  to  come  from  positive  ion  formation,  which  led  to 
NO„  removal  by  both  oxidative  and  reductive  reaction  pathways.  The  relative  humidity  of  the  gas 
played  a  crucial  role  in  the  removal  efficiencies.  The  addition  of  ammonia  led  to  the  formation  of 
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aerosols,  and  a  linear  correlation  was  found  between  N2O  formation  and  the  stoichiometric  ratio  of 
added  ammonia. 

Among  the  disadvantages  of  using  the  electron  beam  process  are  high  cost,  large  physical 
size  of  electron  beam  generators.  X-ray  hazards,  the  necessity  of  ammonia  addition,  and  aerosol 
formation  (Frank  et  al.,  1993). 

Dielectric  Barrier  Discharge 

Dielectric  barrier  discharge  (DBD)  or  silent  discharge  has  been  studied  for  the  removal  of 
SO2  (Sardja  et  al.,  1990;  Dhali  et  al.,  1991)  and  NO,  (Chang  et  al.,  1992).  In  DBD,  a  repetitive 
pulsed  or  AC  voltage  is  applied  in  order  to  sustain  a  discharge  between  two  separated  electrodes 
having  a  dielectric  material  covering  each.  Ions  produced  by  the  discharge  traverse  the  space 
between  the  electrodes,  and  are  then  stored  in  the  surface  of  the  dielectric  materials.  Hence,  a 
sparkless  reactive  plasma  with  a  high  electron  temperature  and  a  low  gas  temperature  is  produced 
which  enhances  the  chemical  reactions  desired  for  SOj/NO,  removal. 

In  DBD,  the  gas  is  exposed  to  many  such  discharge  pulses.  Excitation  of  an  atmospheric 
pressure  plasma  in  this  fashion  is  quite  efficient  because  arcs  between  electrodes  are  prevented  by 
the  short  duration  of  the  individual  current  pulses,  and  the  electron  energy  of  the  plasma  is  well 
matched  with  the  energy  required  to  produce  the  electronic  and  dissociative  states  of  molecular  gases 
(Chang  et  al.,  1992).  Chang  et  al.  (1992)  found  that  under  dry  conditions  the  NO  removal  efficiency 
decreased  with  an  increasing  concentration  of  O2  at  room  temperature.  In  addition,  the  efficiency 
of  NO  removal  increased  in  the  presence  of  water  vapor  in  the  system  at  room  temperature.  At 
elevated  temperatures  (130-160  °C),  increasing  H2O  (>  9%)  concentration  and  decreasing  CO2  (< 
3%)  concentration  resulted  in  increasing  NO  removal  efficiencies.  NO  removal  efficiency  was  found 
to  be  sensitive  to  gas  temperature,  gas  composition,  and  power  deposition. 

Dielectric  discharge  can  also  be  produced  by  utilizing  a  reactor  filled  with  ferroelectric 
material  between  two  metal  electrodes  connected  to  a  high  voltage  AC  power  supply.  An  electric 
field  of  high  intensity  is  formed  around  each  dielectric  pellet,  thereby  producing  high  energy  free 
electrons  (Mizuno  et  al.,  1987). 
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Corona  Discharge 

Corona  discharge  exists  in  several  forms  including  DC  corona,  AC  corona,  or  pulsed  corona, 
of  positive  corona  or  negative  polarity,  depending  on  the  type  and  polarity  of  the  electric  field  and 
on  the  electrode  geometry  (uniform  or  non-uniform).  Pulsed  streamer  corona  discharges  have  been 
considered  for  the  removal  of  pollutants  both  in  the  gas  phase  (Clements  et  al.,  1989,  Dinelli  et  al., 
1990;  Masuda  et  al,  1990)  and  in  the  liquid  phase  (Clements  et  al,  1985;  Sharma  et  al,  1993;  Joshi 
et  al,  1995).  The  use  of  pulsed  corona  for  flue  gas  treatment  is  being  pursued  industrially  in  several 
places  such  as  at  a  coal  burning  thermal  power  station  at  ENEL  (Dinelli  et  al,  1990). 

DC  corona  has  been  used  in  the  removal  of  NO^  in  place  of  an  electron  beam,  and  it  was 
found  that  the  removal  rate  was  enhanced  but  the  power  consumption  was  very  high.  DC  corona 
consumes  high  electric  power  since  it  accelerates  negative  ions  as  well  as  electrons.  Since  oxidizing 
species  such  as  "OH  and  'O  are  produced  in  a  very  short  time  only  by  electrons  with  a  sufficient 
energy,  it  is  imperative  to  accelerate  only  the  electrons,  and  not  the  ions,  in  order  to  reduce  the  power 
consumption.  An  effective  method  is  to  use  a  very  sharp  high  voltage  pulse  with  a  nanosecond  range 
duration  time,  because  this  allows  rapid  acceleration  of  electrons  with  minimum  ion  acceleration 
even  under  ordinary  gas  pressure  (Masuda  et  al,  1990).  Due  to  the  nanosecond  duration  of  the 
applied  pulse,  much  higher  electric  fields  can  be  applied  without  causing  sparking  between  the 
electrodes. 

Simultaneous  removal  of  NO^  and  SO,  from  flue  gas  by  the  application  of  short  duration 
voltage  pulses  was  considered  in  an  electrode  geometry  similar  to  that  of  an  electrostatic  precipitator 
(Dinelli  et  al,  1990).  Nitrogen  oxide  removal  efficiency  of  33%  together  with  a  sulfur  dioxide 
removal  efficiency  of  75%  was  obtained  when  both  contaminants  were  present  at  levels  of  420 
ppmv  and  360  ppmv,  respectively,  in  90  "C  air  containing  ammonia  for  6  Wh/Nra^  energy  input  into 
the  gas.  The  energy  efficiency  of  this  process  is  such  that  with  1  kWh  input  to  the  gas,  a  removal 
of  about  25  g  of  nitrogen  oxide  and  about  70  g  of  sulfur  dioxide  was  achieved  with  the  addition  ot 
40  g  of  ammonia.  The  use  of  hydrated  lime  was  also  considered  as  a  possible  substitute  for  ammonia 
for  the  neutralization  of  the  nitric  and  sulfuric  acid  compounds  to  form  salts.  Optimization  of 
ammonia  injection  and  gas  temperature  along  with  better  electrical  design  of  the  reactors  was 
suggested  to  improve  the  energy  efficiency  . 
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In  a  negative  pulsed  streamer  corona  process,  NO,,  concentration  remained  unchanged  up  tb 
a  critical  peak  field  intensity  for  a  NO-NHj-air  mixture  (NH3 :  180  ppm)  (Masuda  et  al.,  1990). 
Beyond  this  critical  field  intensity,  NO2  was  also  removed,  and  its  removal  depended  upon  the  gas 
residence  time  and  the  amount  of  ammonia  addition.  Ammonia  addition  enhanced  the  removal  of 
NO,  but  not  that  of  NO.  The  efficiency  of  NO  removal  by  using  negative  pulsing  was  found  to  be 
a  function  of  the  specific  power  of  pulsing  P/Q  (corona  power/  gas  flow  rate)  divided  by  the  inverse 
of  the  square  root  of  the  residence  time. 

NO^  removal  from  dry  air  mixed  with  ammonia  in  a  pipe/nozzles-to-plate  electrode  corona 
discharge  system  has  also  been  considered  (Ohkubo  et  al.,  1994).  The  gases  were  introduced  into 
the  discharge  zone  through  the  pipe  electrode  coming  out  of  the  nozzles.  The  addition  of  ammonia 
led  to  the  formation  of  ammonium  nitrate  aerosol  particles  which  were  deposited  on  the  electrodes. 
The  characteristics  of  the  corona  discharge  were  significantly  affected  by  the  addition  of  inert  gases 
such  as  argon.  The  addition  of  these  gases  helped  to  obtain  a  stable  streamer  corona  for  higher  NO„ 
removal  efficiency. 

Mizuno  et  al.  (1995)  noted  that  the  removal  of  NO,  gases  can  be  achieved  by  using  square- 
wave  voltage  pulses  having  frequencies  on  the  order  of  several  hundred  Hz  for  the  energization  of 
cylinder-type  plasma  chemical  reactors.  The  performance  was  noted  to  be  better  at  lower  gas 
temperatures,  and  the  use  of  ethylene  as  an  additive  to  the  flue  gas  significantly  enhanced  the  NO, 
removal  efficiency.  They  observed  the  formation  of  acetic  acid  in  the  analysis  of  the  reaction 
products,  but  it  was  not  quantified. 

Simultaneous  removal  of  NO, ,  SO2  and  fly  ash  (RO,)  using  the  application  of  a  semi-wet 
type  corona  discharge  reactor  has  been  considered  (Okazaki  et  al.,  1995).  The  semi-wet  type 
discharge  reactor  consisted  of  a  cylindrical  pyrex  glass  tube  and  a  stainless  steel  discharge  wire 
electrode  at  the  center,  with  a  water  film  covering  the  inner  surface  of  the  glass  tube.  For  an  initial 
concentration  of  300  ppm  of  NO  and  500  ppm  of  SO2,  90%  of  NO,  and  nearly  100%  of  SO2  were 
removed  at  the  same  time,  even  in  a  gas  containing  a  high  oxygen  concentration  of  about  15%. 

Tas  (1995)  investigated  the  plasma-induced  conversion  of  NO  as  a  function  of  several 
parameters  including  the  composition  of  the  gas  mixture,  the  presence  of  catalyst  materials,  and  the 
type  of  catalyst  material.  The  energy  efficiency  of  the  conversion  of  NO  to  NO2  can  be  increased 
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by  a  factor  of  2  for  pure  bulk  gas  of  helium  or  nitrogen  with  1000  ppm  of  NO  by  using  a  silica 
catalyst.  The  energy  input  was  varied  by  controlling  the  pulsing  frequency  of  the  discharge  in  the 
reactor.  The  addition  of  oxygen  to  the  bulk  gas  depressed  the  dissociative  conversion  of  NO,  and 
at  a  high  power  input  and  for  an  initial  concentration  of  1000  ppm,  equation  (1)  corresponded  well 
with  the  measured  conversion  of  NO.  With  the  addition  of  water,  NO2  was  expected  to  be 
converted  to  HNO3,  however  HNO3  was  not  detected. 

Figure  1  compares  the  energy  efficiency  of  NO  removal  at  high  concentrations  of  NO  by 
Mizuno  et  al.  (1995)  and  Tas  (1995),  and  Figure  2  shows  a  summary  of  energy  efficiency  of  NO 
removal  at  low  concentrations  of  NO  by  Masuda  et  al.  (1990)  and  Ohkubho  et  al.  (1994).  As  seen 
in  Figure  1,  the  addition  of  500  ppm  of  ethylene  increased  the  NO  removal  for  the  same  power  per 
unit  volume  divided  by  the  square  root  of  time  in  dry  air.  The  data  of  Mizuno  et  al.  ( 1 995)  and  Tas 
(1995)  for  the  case  of  NO  removal  in  the  presence  of  water  are  very  close. 


1400.0 

1200-0 

1000.0 

>• 

o 

e 

^  800.0 
O 

£  600.0 
o. 

O- 

400.0 

200.0 

0.0 


■  Mizuno  (10%  O2) 

•  Mizuno  (10%  O2,  SOOppm 

♦  Mizuno  (10%  O2.  4-. 5  or  9%  HgO) 
▼  Tas(5%  O2,  1 1  °/o  H2O  ) 


200.0 


0.0  100.0 

(Power/Volumo).  Tlme  ’^  ((WH/m  ).s©c‘  ) 


300.0 


Figure  1.  Summary  of  Energy  Efficiency  of  NO  Removal  at  High  Concentrations. 


Figure  2  shows  that  the  NO  removal  in  the  presence  of  ammonia  considered  by  Masuda  and 
Nakao  (1990)  was  far  different  from  that  observed  by  Ohkubho  et  al  (1994).  In  all  of  the  above 
cases  where  NO  removal  was  considered  by  pulsed  streamer  corona  discharge,  there  were  certain 
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dissimilarities  such  as  the  nature  of  the  electrode  and  the  polarity  of  the  field.  No  theoretical 
justification  for  this  type  of  plot  is  available  at  present,  however  it  serves  as  a  convenient  means  to 

compare  data  under  a  variety  of  power  conditions. 

Modeling  of  the  chemical  reactions  taking  place  during  the  pulsed  streamer  corona  process 
has  been  developed  in  order  to  understand  the  evolution  of  the  concentrations  of  different  chemical 
species  (Penetrante  et  al.,  1993).  For  a  power  input  of  0.05  W/cm^  and  for  an  initial  concentration 
of  400  ppm  NO  at  300  K,  200  ppm  of  the  initial  NO  was  predicted  to  be  converted  back  to  N,  and 
200  ppm  to  NO2.  Of  this  200  ppm  of  NO  2, 1 80  ppm  was  converted  to  N  ,0  and  20  ppm  to  N  2O  5 
When  the  temperature  was  increased  to  400  K,  only  100  ppm  of  the  initial  NO  was  converted  to  NO, 
and  another  100  ppm  was  converted  to  N2O.  With  the  addition  of  10%  water  and  15%  CO2,  from 
an  initial  400  ppm  NO  about  120  ppm  was  converted  to  HNO3  and  50  ppm  was  converted  to  N2O. 
The  ozone  produced  was  about  50  ppm  and  the  residual  NO^  present  was  around  10  ppm.  This 
model  prediction,  however,  has  not  yet  been  confirmed  by  experimental  results. 
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Figure  2.  Summary  of  Energy  Efficiency  for  NO  Removal  at  Low  Concentration. 


Recent  experimental  studies  by  Wolf  et  al.  (1996)  of  dielectric  barrier  discharges  have  shown 
that  NO  removal  depends  not  only  on  the  power  density  delivered  to  the  reactor,  but  also  on  the 
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absolute  value  of  the  applied  voltage.  This  agrees  with  the  work  of  Masuda  et  al.  (1990)  and 
Mizunoetal.  (1995)  who  found  the  same  result  for  pulsed  corona  reactors.  Penetrante  et  al.  (1996) 
have  claimed  that  the  performance  of  dielectric  barrier  discharge  reactors  should  be  the  same  as  that 
of  pulsed  corona  reactors,  and  that  the  performance  of  each  should  only  depend  upon  the  power 
density  delivered  to  the  reactor  and  not  upon  the  field  strength  or  characteristic  shape  of  the  voltage 
pulse.  They  showed  that  NO  removal  versus  power  density  was  the  same  for  specially  designed 
corona  and  dielectric  barrier  reactors.  Wolf  et  al.  (1996)  constructed  a  dielectric  barrier  reactor  that 
would  operate  at  higher  fields  and  found  that  the  power  density  was  lowered  for  a  given  NO 
removal. 

Physical  Characteristics  of  the  Pulsed  Corona  Discharge 

Different  types  of  corona  discharges  can  occur  in  non-uniform  electric  fields,  depending  upon 
the  polarity  and  the  namre  of  the  applied  electric  field.  Positive  pulsed  corona  proceeds  from  corona 
onset  to  a  burst  pulse  corona  to  streamer  corona,  then  to  glow  corona  and,  finally,  spark  discharge 
as  the  applied  voltage  is  increased.  In  the  case  of  negative  pulsed  corona,  the  initial  form  is  Trichel 
pulse  corona,  followed  by  pulseless  corona  and  then  spark  discharge  as  the  voltage  is  increased.  A 
Trichel  pulse  refers  to  the  regular  current  pulse  produced  when  the  voltage  applied  to  a  negative 
electrode  is  gradually  increased  to  produce  ionization.  These  pulses  are  extremely  regular  in  both 
magnitude  and  repetition  rate.  The  frequency  of  the  Trichel  pulse  repetition  starts  with  2  kHz  near 
onset,  and  it  may  attain  a  few  MHz  before  the  onset  of  glow  corona  (Naseer,  1971).  Positive  corona 
depends  more  on  photoionization  for  its  propagation,  while  the  negative  corona  propagates  by  impact 
ionization  of  the  gas  molecules  (Chang  et  al.,  1993). 

In  pulsed  streamer  corona,  the  ionization  grows  from  the  high  voltage  electrode  into  a 
background  field  which  is  purely  geometric.  The  discharge  assumes  the  form  of  a  number  of 
branched,  tree-like  streamers,  whose  starting  points  are  discrete  and  distributed  over  the  highly 
stressed  region  with  a  specific  density  (Dinelli  et  al,  1990).  Photoionization  is  generally  assumed 
to  be  the  mechanism  that  supplies  secondary  seed  electrons  and  triggers  the  avalanches. 

The  streamer  formation  in  a  DC  anode  (positive)  corona  can  be  described  as  a  Townsend 
avalanche  which  is  a  conceptual  description  of  the  electron  multiplication  processes.  The  process 
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begins  with  free  electrons  in  a  strong  electric  field  gaining  energy  to  collide  with  gas  molecules  and 
lose  energy  by  a  variety  of  collision  process  (e.g.  elastic  scattering,  excitation,  ionization),  thereby 
producing  secondary  electrons  and  positive  ions  (Creyghton,  1 994).  The  electrons  produced  through 
collision  further  produce  electrons  through  secondary  reactions,  resulting  in  an  electron  avalanche. 
As  a  result,  the  head  of  the  electron  avalanche  and  the  positive  ions,  which  are  left  behind  the 
avalanche,  form  a  sort  of  dipole.  The  electron  avalanche  is  said  to  be  "critical"  when  the  dipole  field 
reaches  a  value  comparable  to  the  external  applied  field.  The  electric  field  is  enhanced  in  front  of 
the  head  of  the  electron  avalanche,  while  the  electric  field  is  reduced  in  between  the  avalanche  head 
and  the  positive  ions  which  are  left  behind.  Still  further  back  in  the  avalanche,  the  electric  field 

between  the  positive  ions  and  the  cathode  is  enhanced. 

The  electron  avalanche  has  an  increased  velocity  and  light  emission  intensity.  These 
avalanches  are  referred  to  as  anode  directed  streamers  (ADS).  Also,  as  the  electrons  propagate 
towards  the  anode  there  is  a  buildup  of  positive  ions.  In  the  presence  of  a  space  charge  field  and 
initial  free  electrons  at  the  cathode  side  of  the  avalanche,  cathode  directed  streamers  (CDS)  can  be 
formed.  The  CDS  propagate  as  a  result  of  secondary  avalanches  which  neutralize  the  positive 
concentration  of  the  previous  avalanches,  leaving  a  new  positive  ion  concentration  behind  at  a 
shorter  distance  from  the  cathode  (Creyghton,  1994).  During  the  primary  avalanche  propagation 
excited  gas  molecules  are  also  produced  which  act  as  a  source  for  electrons.  These  electrons  serve 
as  a  source  for  secondary  avalanches  required  for  CDS  propagation. 

Streamer  discharges  are  hence  highly  localized  "space  charge  waves"  which  enhance  the 
applied  field  in  front  of  the  wave  (active  region)  and  propagate  because  of  electron  avalanche 
formation  in  this  high  field.  Along  the  track  of  the  wave  there  remains  a  weakly  ionized  plasma 
filament  (passive  region),  along  which  the  conduction  current  flows  to  the  high  voltage  electrode, 
supplying  the  energy  for  streamer  advancement  (Nasser,  1971).  In  the  case  of  corona  at  impulse 
voltages,  the  ionic  space  charge  drift  and  accumulation  can  be  neglected  due  to  the  short  duration 
of  the  pulse  width.  In  the  case  of  non-uniform  electrode  geometry,  the  streamer  growth  assumes 
radial  directions,  increasing  the  number  of  branches. 

The  major  characteristics  of  streamers  include  the  length,  number,  speed  of  propagation, 
width,  number  of  electrons  and  ions,  time  lag,  and  electric  field  distribution.  Typical  values  of 
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streamer  width  are  in  the  range  of  25  to  100  ^tm  and  the  speed  of  propagation  is  in  the  range  of  0.7- 
0.9  mm/ns  (Creyghton,  1994).  Recent  experimental  and  theoretical  studies  have  also  shown  that  the 
propagation  characteristics  of  positive  streamers  depend  very  much  upon  the  degree  of 
electronegativity  of  the  carrier  gas.  In  weakly  electronegative  gases  like  atmospheric  air,  the 
streamers  appear  much  longer  and  more  branched  than  in  strongly  electronegative  gas  such  as  SF5. 
Hence,  with  large  concentrations  of  CO2  and  H2O,  the  discharge  characteristics  are  very  near  to  those 
of  highly  electronegative  gases  (Gallimberti,  1988). 

Modeling  Pulsed  Streamer  Corona  Discharge 

In  order  to  model  the  streamer  discharge  in  a  pulsed  corona,  an  understanding  of  the 
elementary  processes  such  as  ionization,  attachment,  and  photoionization  is  essential. 

A  1-D  streamer  propagation  model  has  been  developed  in  which  the  streamer  is  assumed  to 
possess  cylindrical  symmetry  (Gallimberti,  1972, 1988).  The  densities  of  the  charged  particles  along 
the  central  axis  (x-axis)  of  the  streamer  can  be  calculated  solving  the  continuity  equations,  i.e.,  for 
electrons,  positive  ions,  and  negative  ions.  It  is  possible  to  solve  for  the  electric  field  by  solving  the 
Poisson  equation  once  the  densities  of  the  charged  particles  are  known; 


where  E  is  the  electric  field,  np  is  the  concentration  of  positive  ions,  n^  is  the  concentration  of 
electrons,  n„  is  the  concentration  of  negative  ions  and  is  the  permittivity  of  the  gas  medium 
considered. 

2-D  models,  in  which  the  radial  distribution  of  charged  particles  is  calculated,  have  been 
performed  to  simulate  streamer  formation  in  a  uniform  electric  field  (Wang  et  al.,  1990).  A  1-D 
streamer  model  with  an  approximation  of  the  Poisson  equation  has  been  used  to  calculate  the 
streamer  properties  (Creyghton,  1994).  The  current  pulse,  net  charge,  and  energy  loss  have  also 
been  calculated  for  an  impulse  corona  (Gallimberti,  1988). 

The  solution  of  the  Boltzmann  equation  is  necessary  to  calculate  the  electron  transport 
coefficients  and  kinetic  rate  coefficients  for  various  inelastic  electron-gas  collision  reactions.  The 
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Boltzmann  equation  is  the  continuity  equation  for  electrons  in  six-dimensional  phase  space  and  it 
describes  the  time  evolution  of  the  electron  energy  distribution  function.  The  most  common 
technique  used  to  solve  the  Boltzmann  equation  has  been  the  moment  method  in  which  a  spherical 
harmonic  expansion  (or  Legendre  series)  of  the  angular  dependence  of  the  velocity  is  assumed  with 
truncation  after  the  first  two  terms  (Huxley  at  al.,  1974).  The  electron  transport  properties  were 
significantly  affected  by  the  composition  of  gas  considered. 

The  kinetic  rate  coefficients  for  the  electron-gas  reactions  are  functions  of  the  electron  energy 
distribution  function.  The  kinetic  rate  coefficients  can  be  calculated  by 

0  ^ 


where  k.  is  the  kinetic  rate  constant  for  the  electron-gas  reaction,  s  is  the  electron  energy,  m,  is  the 

ClJ 

mass  of  an  electron,  Gy  is  the  collisonal  cross  section,  and  f(s)  is  the  electron  energy  distribution 
function. 

The  collision  cross  section  is  defined  as  a  measure  of  the  probability  of  a  given  type  of 
interaction  between  a  gas  molecule  and  an  electron.  The  cross  section  is  a  function  of  the  relative 
velocity  of  particles  and  the  energy  of  the  electrons.  It  varies  over  several  orders  of  magnitude 
according  to  the  type  of  interaction  being  considered. 

The  densities  of  chemical  reactive  dissociation  products  have  been  calculated  using  a 
numerical  model  for  streamer-induced  chemical  kinetics  as  discussed  above  (Alekseev  et  al.,  1993, 
McFarlane  et  al.,  1991 ;  Maroude  et  al.,  1996).  Modeling  of  the  electrostatic  corona  discharge  reactor 
has  been  done  by  solving  the  Boltzmann  equation  for  the  electron  energy  distribution  function,  which 
in  turn  was  used  to  determine  the  kinetic  rate  constants  for  the  electron  gas  collisional  reactions 
(Mukkavilli  et  al.,  1988).  For  a  dry  air  system,  58  reactions  were  considered,  and  96  reactions  were 
considered  for  the  wet  air  system.  Increase  of  moisture  in  the  feed  air  had  a  negligible  effect  on  the 
input  gases  (N,  and  O2),  but  it  was  found  to  reduce  the  ozone  concentration  by  more  than  two  orders 
of  magnitude  and  that  of  NO^  by  a  factor  of  two. 
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Ozone  generation  in  pulsed  corona  reactors  having  point-to-plane  and  wire-cylinder 
geometries  have  been  considered  with  discrete  injection  of  electrons  and  also  with  time-varying 
electronic  density  (Pignolet  et  al,  1990;  Loiseau  et  al.,  1994).  The  shape  of  the  electronic  injection 
time  profile  did  not  affect  the  concentration  profile  of  ozone  production.  The  temperature  gradient 
between  the  electrodes  in  the  wire-to-cylinder  geometry  had  a  drastic  effect  on  the  electron  densities 
near  the  wire  directly  through  ozone  destruction  at  high  temperatures  and  indirectly  by  viscosity, 
which  affects  the  velocity  profile  (Loiseau  et  al.,  1994). 

Past  work  in  the  field  of  pulsed  streamer  corona  has  focused  on  either  primarily  experimental 
results  or  on  purely  theoretical  modeling.  The  present  work  will  try  to  combine  both  experimental 
and  computational  aspects  in  order  to  develop  an  understanding  of  the  basic  chemistry  taking  place 
during  the  pulsed  streamer  corona  process.  Removal  of  NO  from  gases  of  various  compositions 
including  pure  nitrogen  and  air,  and  the  addition  of  water  vapor,  ethylene  and  carbon  monoxide  will 
be  considered.  A  kinetic  model  will  be  utilized  in  order  to  develop  a  quantitative  understanding 
of  the  reaction  pathways  occurring  in  the  system. 
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SECTION  II 

EXPERIMENTAL  APPARATUS  AND  PROCEDURES 


Non-thermal  plasma  reactors  operated  at  low  temperature,  low  power  consumption,  and 
atmospheric  pressure  can  be  classified  by  their  physical  geometry  and  by  their  method  of  electrical 
energization.  Among  the  various  types  are  point-to-plane,  point-to-point,  packed  bed,  coaxial,  and 
wire-cylinder  electrode  geometries.  The  methods  of  energization  include  direct  current,  alternating 
current,  or  pulsed  current.  In  this  study,  pulsed  power  was  delivered  to  a  wire-cylinder  geometry 
reactor. 

A.  POWER  SUPPLY 

The  pulsed  power  supply  used  in  this  investigation  is  shown  in  Figure  3.  It  consists  of  two 
major  units:  1)  a  high  voltage  transformer/rectifier  (T-R)  set  and  2)  a  pulse  generating  apparatus. 
The  power  supply  used  was  similar  to  the  configuration  used  by  Clements  et  al.  (1989).  The  high 
voltage  DC  T-R  set  was  modified  by  eliminating  the  internal  rectifier,  thereby  producing  0-100  kV 
AC  at  0-28  mA.  Output  from  the  T-R  set  was  connected  to  the  pulse  generation  circuit,  where  the 
current  passes  through  a  high  voltage  resistor  (333  kQ)  and  then  through  a  diode  array  which  acts 
as  a  half-wave  rectifier.  The  output  from  the  rectifier  charged  a  bank  of  capacitors  (2700  pF)  during 
the  charging  cycle.  A  mechanical  rotating  spark  gap  having  two  fixed  stainless  steel  spherical 
electrodes  and  a  rotating  rod  electrode  was  employed  to  discharge  the  capacitor  bank  twice  per  shaft 
revolution.  s>Tichronized  with  the  AC  current  (60  Hz).  From  the  output  of  the  spark  gap  the  circuit 
had  two  parallel  pathways.  One  of  them  passed  through  a  load  resistor  (300  kQ)  which  enabled  the 
gap  to  fire.  The  other  was  connected  to  the  central  wire  within  the  cylindrical  reactor,  from  which 
the  discharge  emanated.  The  generation  of  pulsed  corona  discharge  by  a  rotating  spark  gap  was  a 
source  of  electromagnetic  interference,  therefore  the  pulse  generating  apparatus  was  encased  in  an 
aluminum  box  to  reduce  RF  noise. 


16 


HVAC  =  0  -  100  kV,  0-28  mA 
=  333  kO 
Rg  =  300  0 

=  60  kV,  12  A  max 

=  2700  pF 

=  Rotating  Spark  Gap,  1800  RPM,  60  discharges  /  sec 
=  Strobe  Lamp 


Figure  3  ;  Pulsed  Power  Supply  Diagram 
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B.  REACTOR  SPECIFICATIONS 

The  pulsed  streamer  corona  reactor  used  in  this  study  was  a  stainless  steel,  4  inches  in 
diameter  and  18  inches  long,  cylindrical  tube.  The  high  voltage  electrode  (0.109  inches  diameter, 
12  inches  long)  was  suspended  concentric  to  the  grounded  outer  cylinder  resulting  in  an  electrode 
separation  of  2  inches.  The  active  pulsed  corona  treatment  volume  of  the  reactor  was  2.45  /.  The 
typical  gas  flow  rate,  gas  velocity,  and  residence  time  in  the  reactor  are  from  2.5  to  40  //min,  from 
321  to  5134  cm/s,  and  from  3.6  to  60.0  seconds,  respectively.  The  high  voltage  ceramic  feed¬ 
through  conductors  as  well  as  the  center  electrode  were  fabricated  of  #  304  stainless  steel  to  prevent 
any  corrosion  and  rust  formation  inside  the  reactor. 

C.  ANALYTICAL  INSTRUMENTATION 

Analysis  of  the  outlet  gases  from  the  reactor  included  the  measurement  of  NO,  NO..  O-,, 
N^O,  and  hydrocarbon  breakdown  products.  Nitrogen  oxides  (NO  and  NO2)  were  measured  using 
a  Thermo-Environmental  Instruments,  Inc.,  Model  42H  chemiluminescence  NO,;  analyzer.  Ozone 
was  measured  using  the  KI  reaction  test.  N^O  and  hydrocarbon  breakdown  products  were  measured 
using  a  Hewlett-Packard  GC/MS  instrument  located  in  the  Department  of  Chemistry  at  Florida  State 
University.  HPLC  analysis  was  also  performed  on  samples  obtained  by  sparging  the  reactor  exhaust 
gas  through  a  flask  of  HPLC  carrier  solvent  (i.e.,  phosphoric  acid  buffer  in  de-ionized  water).  HPLC 
analysis  was  used  to  determine  the  breakdown  products  of  the  ethylene  when  NO  removal  was 
considered  in  the  presence  of  ethylene  in  humid  air.  The  pulse  waveform  characteristics  were 
measured  using  a  Tektronix  TDX  460  four  channel  digitizing  oscilloscope  with  a  P6015A  lOOOx  3.0 
pF,  1000  MO  voltage  probe  and  a  current  probe  (P6016).  A  Cole-Parmer  TriSense 
velocity/temperature/humidity  measurement  system  was  used  to  measure  the  gas  velociw. 
temperature,  and  humidity  in  the  gas  system.  The  analytical  instrumentation,  the  reactor  setup,  and 
the  gas  feed  system  are  shown  in  Figure  4a. 

D.  FEED  GAS  SYSTEM 

The  feed  gas  system  (see  Figure  4a)  consisted  of:  1)  a  compressed  air  source  and  pine 
nitrogen  gas  cylinders  which  act  as  carrier  gases,  2)  gas  cylinders  and  high  pressure  regulators 
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containing  1%  NO  in  dry  Nj,  1%  NOj  in  dry  Nj,  pure  C2H4,  and  pure  CO,  3)  gas  cylinders  and  high 
pressure  regulators  containing  calibration  gases  of  100  ppm  NO  in  Nj  and  113  ppm  NO2  in  N2, 4) 
a  water  tank  through  which  the  air  and  the  other  gas  can  be  bubbled  for  humidification,  5)  flow 
meters,  valves,  and  mass  flow  controllers,  and  6)  a  stainless  steel  gas  mixing  chamber  upstream  from 
the  reactor. 

The  gas  cylinders  and  the  regulators  were  manufactured  by  Air  Products  Inc.  The  water  tank 
used  to  humidify  the  air  was  26  inches  long  and  7  inches  in  diameter.  The  air  was  introduced  at  the 
bottom  of  the  tank,  and  was  made  to  bubble  through  the  water  eolumn.  Flow  meters  used  to  monitor 
the  flow  rate  of  the  air  into  the  reactor  were  purchased  from  Dwyer,  Inc.  Mass  flow  controllers  used 
to  control  the  flow  of  air  and  the  trace  gases  sent  into  the  reactor  were  purchased  from  MKS 
Instruments  Inc.  Two  controllers  of  Type  1259C  were  used  to  control  the  flow  of  NO  and  ethylene 
gas.  A  mass  flow  controller  of  Type  1559A  was  used  to  control  the  amount  of  air  or  N2  flowing  into 
the  reactor.  Trace  gases  such  as  NO,  ethylene,  and  CO  were  mixed  with  the  air  in  a  stainless  steel 
mixing  chamber  before  passing  into  the  reactor.  The  mixing  chamber  was  23.5  in.  long  and  had  a 
dijimeter  of  2.5  inches. 

E.  EXPERIMENTAL  PROCEDURES 

In  order  to  conduct  a  typical  experiment  to  consider  NO  removal  at  a  particular  voltage  and 
residenee  time,  the  NO^  monitor  as  well  as  the  mass  flow  controllers  were  turned  on  at  least  2  hours 
before  the  start  of  the  experiment  to  allow  them  to  stabilize.  Calibration  of  the  NO^  monitor  was 
performed  by  passing  the  calibration  gases  of  known  concentrations  of  NO  and  NO2  into  the  monitor 
and  calibrating  the  monitor  reading  to  the  corresponding  value.  Once  the  mass  flow  controllers 
stabilize,  the  readings  were  zeroed,  if  neeessary. 

The  feed  air  used  in  the  experiments  was  produced  by  a  central  compressed  air  source  from 
the  College  of  Engineering  building.  A  dedicated  regulator/valve  filter  assembly  was  adjusted  so 
that  the  pressure  of  the  air  delivered  to  the  gas  distribution  system  was  maintained  at  20  psig.  The 
air  flow  rate  through  the  reactor  was  then  adjusted  with  needle  valves  and  with  the  mass  flow 
controller  corresponding  to  the  required  gas  residence  time  in  the  reactor.  The  two  spherical 
electrodes  and  the  rod  electrode  on  the  rotating  spark  gap  tended  to  corrode  during  the  course  of  an 


experiment,  so  they  were  cleaned  using  ethanol  and  sandpaper  before  running  every  experiment. 
The  rotating  mechanical  spark  gap  was  engaged  and  adjusted  so  that  it  was  aligned  perpendicular 
to  the  discharge  position  when  the  strobe  lamp  was  in  line  frequency  mode  (60  Hz).  This  ensured 
full  charging  of  the  capacitors  by  the  power  supply. 

The  feed  air  delivered  to  the  reactor  was  first  passed  through  a  drier  (calcium  sulfate  packed 
columns)  to  remove  any  residual  water  vapor  present  in  the  compressed  air.  In  cases  where  water 
vapor  was  required  for  an  experiment  in  the  system,  air  was  bubbled  through  a  water  tank.  The  flow 
rate  of  all  feed  and  trace  gases  were  continually  monitored  by  mass  flow  controllers  and/or  by  flow 
meters. 

To  measure  the  voltage  waveform  of  the  pulsed  voltage,  the  high  voltage  probe  was 
connected  to  the  central  wire  electrode  to  which  the  high  voltage  input  was  connected.  The  current 
probe  was  attached  across  a  wire  connecting  the  outer  cylinder  of  the  reactor  to  the  ground.  The 
exact  locations  at  which  the  voltage  and  the  current  were  measured  are  shown  in  Figure  4.  Both  the 
probes  were  connected  to  the  oscilloscope  to  measure  the  different  characteristics  of  the  waveforms. 
Once  the  flow  rate  of  all  the  gases  had  stabilized,  the  pulsed  corona  discharge  was  applied  to  the 
reactor. 

The  pulsed  voltage,  current,  and  power  waveform  characteristics  were  noted  from  the 
oscilloscope,  and  adequate  time  was  allowed  for  all  the  analytical  instruments  to  stabilize.  The 
exhaust  gases  coming  out  of  the  reactor  were  passed  through  the  NO,;  monitor  to  record  the 
concentration  of  NO  and  NOj.  For  trace  gas  sample  analysis,  the  outlet  gases  from  the  reactor  were 
passed  through  a  glass  sampling  vessel.  The  gases  in  the  sampling  vessel  were  analyzed  by  GC/MS 
for  the  breakdown  products  of  NO,;  gases  and  also  for  the  concentration  of  ethylene  in  experiments 
for  which  the  ethylene  was  present  in  the  feed  gas. 

For  liquid  sample  analysis,  the  exhaust  gases  were  sparged  through  water  in  a  flask.  These 
water  samples  were  then  analyzed  by  HPLC  and  ion  chromatography.  HPLC  analysis  was  performed 
on  a  Perkin  Elmer  HPLC  unit  using  a  Supelcogel  C-610H  organic  acid  column  30  cm  by  7.8  mm. 
The  flow  rate  was  0.5  ml/min  with  a  mobile  phase  of  0.1%  H3PO4.  A  UV  detector  set  at  210  nm  was 
used  and  samples  were  injected  through  a  300  pi  sample  loop.  HPLC  Grade  Optima  water  was  used 
for  the  mobile  phase  preparation  and  for  the  collection  of  the  samples  from  the  corona  reactor. 
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Exhaust  gas  from  the  corona  reactor  was  bubbled  through  approximately  200  ml  of  the  phosphoric 
acid  mobile  phase  for  specified  time  intervals  of  20  to  30  minutes.  Samples  of  this  fluid  were 
analyzed  by  HPLC.  The  exhaust  gases  used  for  this  analysis  were  obtained  from  the  corona  with 
operating  conditions:  1)  50  kV,  2)  residence  time  44  seconds,  3)  100  ppm  NO,  4)  500  ppm  C2H4,  5) 
carrier  gas  of  compressed  air  humidified  to  85.5%  relative  humidity  at  75  °C. 

Ozone  concentration  was  determined  by  using  the  KI  method  specified  in  Greenberg  et  al. 
(1992).  This  involved  the  use  of  two  traps  in  series  each  containing  200  mL  of  2%  KI  solution  at 
the  output  of  the  corona  reactor.  The  entire  output  flow  from  the  reactor  was  sparged  into  each  trap 
in  series.  The  traps  were  exposed  to  the  output  stream  for  10  minutes.  After  exposure  both  solutions 
were  transferred  to  separate  beakers  and  titrated  against  0.05  N  Na2S203  solution.  When  the 
solutions  lose  their  yellow  color,  1  to  2  mL  of  starch  indicator  solution  was  added.  Titration  was 
continued  until  the  solution  was  clear.  Ozone  concentration  is  calculated  using  the  equations  given 
by  Greenberg  et  al.  (1992). 

Once  the  readings  of  the  1)  NO^  concentration,  2)  rise  time,  3)  peak  voltage,  4)  pulse  width 
of  the  voltage  waveform,  5)  peak  current,  and  6)  pulse  width  of  the  current  waveform  were 
measured,  the  power  to  the  reactor  was  turned  off  and  the  reactor  and  the  spark  gap  were  grounded 
with  a  grounding  rod.  The  valves  of  the  gas  cylinders  were  closed  and  the  reactor  was  then  purged 
for  several  minutes  with  clean,  dry  compressed  air. 

F.  LIQUID  PHASE  REACTOR  SPECIFICATIONS 

A  new  liquid  phase  reactor  was  designed  and  constructed  following  the  general  design  of 
reactors  used  in  previous  studies  (Sharma  et  al.,  1993,  Joshi  et  al.,  1995).  The  reactor  consists  of  a 
I  liter  jacketed  glass  vessel  shown  in  Figure  4b.  A  glass  u-tube  runs  from  outside  the  jacketed 
reactor  to  a  male  ground  glass  connector  on  the  inside  of  the  reactor.  A  high  voltage  wire  will  be 
introduced  into  the  reactor  through  this  tube;  the  wire  will  connect  with  a  stainless  steel  discharge 
needle  facing  upwards.  Opposite  the  needle,  a  pancake-shaped  grounded  electrode  is  suspended 
from  a  ground  glass  fitting  centered  in  the  top  of  the  reactor  lid.  Another  ground  glass  fitting  will 
allow  gas  to  be  introduced  into  the  reactor,  while  a  third  fitting  will  be  impermeable. 

The  improvements  in  the  reactor  design  from  the  previous  reactors  used  in  our  laboratory  are 
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Figure  4b.  Liquid  Phase  Reactor  and  Electrode  Tip  Design. 


23 


the  following.  Reactor  volume  has  been  increased  from  550  ml  to  one  liter,  and  glass  is  used  as  the 
vessel  material.  The  high  voltage  input  wire  will  be  introduced  from  the  side  of  the  reactor  and  will 
not  be  in  close  proximity  to  the  ground  electrode.  This  will  also  allow  for  a  perfectly  round  ground 
electrode  that  will  not  have  any  sharp  edges.  The  addition  of  the  water  jacket  will  allow  for  easier 
control  of  the  reactor  temperature.  Finally,  the  new  design  of  the  electrode  tip  will  allow  for  easy 
removal  of  the  electrode  from  the  reactor  and  replacement  of  the  discharge  electrode. 
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SECTION  III 
REACTION  MODELING 


In  a  pulsed  corona  discharge  in  a  gaseous  medium,  free  electrons  gain  energy  from  an 
imposed  electric  field  and  lose  energy  through  collisions  with  neutral  gas  molecules.  This 
transfer  of  energy  to  the  gas  molecules  leads  to  the  formation  of  a  variety  of  new  species 
including  ions,  metastables,  atoms,  and  free  radicals.  These  products  are  chemically  active  and 
react  with  the  other  gas  molecules  to  form  new  stable  compounds.  The  modeling  of  the  reactor  is 
a  complicated  process  and  certain  assumptions  are  made  to  analyze  the  system. 

The  kinetic  rate  constants  for  the  electron-gas  reactions  depend  upon  the  imposed  electric 
field  and  also  on  the  gas  composition.  In  this  work,  a  model  was  constructed  with  the  help  of 
experimental  data  to  consider  the  concentration  of  different  molecular  species  during  the  pulsed 
corona  discharge.  A  global  approach  was  adopted  where  only  dissociation  reactions  of  major 
constituents  of  the  gas  were  considered  in  the  reaction  model.  The  kinetic  rate  constants  for  the 
dissociation  reactions  were  obtained  for  the  experimental  data  on  NO  removal  and  NO, 
formation  and  removal  in  various  gas  streams. 

The  modeling  effort  here  accounted  for  the  chemistry  in  pulsed  streamer  corona  discharge 
by  considering  two  chemical  reaction  sets.  During  the  time  frame  of  pulsing  (pulse-on  period), 
the  electron  concentration  was  assumed  to  be  fixed,  and  the  reaction  set  included  the  gas 
molecule  dissociation  reactions  in  order  to  calculate  the  concentration  profile  of  different  species. 
During  the  time  between  the  pulsing  (pulse-off  period),  the  electron  concentration  was  assumed 
to  be  zero,  and  so  the  reaction  set  which  did  not  have  the  dissociation  reactions  was  used  to 
calculate  the  concentration  profile.  These  two  reaction  sets  were  interchanged  repeatedly  to 
simulate  the  pulse-on,  pulse-off  aspect  of  the  discharge. 

The  pulse  width  of  the  voltage  pulse,  which  was  the  time  at  the  half  maximum,  varied 
between  500-700  nanoseconds.  The  total  pulse  duration  from  the  inception  of  the  rise  time  to  the 
complete  decay  of  the  voltage  peak  was  1  microsecond.  In  the  present  modeling,  the  pulse  width 
was  assumed  to  be  1  microsecond  and  the  time  between  the  two  pulses  was  assumed  to  be  16.667 
milliseconds,  corresponding  to  a  pulse  frequency  of  60  Hz.  Hence  the  reaction  set  Avith  electron 
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dissociation  reactions  was  used  to  calculate  the  concentration  profile  for  1  microsecond  followed 
by  the  reaction  set  which  did  not  include  the  dissociation  reactions  to  calculate  the  concentration 
profiles  for  the  next  16.667  milliseconds  repeatedly  until  the  final  total  time  equaled  the 
residence  time  in  the  reactor. 

This  method  was  similar  to  the  model  used  by  Loiseau  et  al.  (1994)  for  computing  ozone 
production  in  a  cylindrical  oxygen-fed  ozonizer.  For  example,  Figure  5  shows  the  concentration 
profile  of  ozone  computed  by  the  code  used  in  the  present  study  that  alternately  employed  the 
two  reaction  sets  to  simulate  the  pulse-on  and  pulse  off  aspect  of  the  discharge.  As  in  the  figure, 
during  the  pulse-on  period  the  oxygen  radicals  are  produced  and  hence  ozone  formation 
increased  rapidly,  and  during  the  pulse-off  period  the  ozone  production  level  was  flat. 


Figure  5.  Ozone  Concentration  Profile  to  Illustrate  Modeling  Method. 

The  following  were  the  assumptions  made  for  the  model  used  here: 

1 .  Isothermal  system.  There  is  no  appreciable  heating  of  the  gas  during  the  residence  time  in  the 
reactor. 

2.  Plug-flow  model.  The  model  does  not  account  for  the  axial  and  radial  dispersion  in  the 
reactor. 

3.  There  is  no  variation  in  the  concentration  of  electrons  during  the  pulse-on  period  either  in 
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space  or  in  time. 

4.  Reactions  involving  metastables,  ions,  and  photo-ionization  are  neglected. 

5.  The  model  does  not  account  for  the  variation  in  electric  field  strength  from  the  anode  wire  to 
the  cathode  cylinder. 

The  molar  species  continuity  equations  for  all  the  species  is  given  by: 


dc 

dt 


‘-*u-Vc,-Dfc,*Y.  "A 


(4) 


where  c,  is  the  concentration  of  species  i,  t  is  the  time,  u  is  the  molar  average  velocity  vector,  D|  is 
the  diffusion  coefficient,  Vy  is  the  stoichiometric  coefficient  of  the  ith  species  in  the  jth  reaction  and 
i?y.  is  the  reaction  rate.  For  a  plug  flow  reactor  model  without  axial  dispersion  or  molecular  diffusion 
the  equation  reduces  to: 


w 


(5) 


In  terms  of  residence  time  (t=x/uJ  in  the  reactor,  the  equation  becomes: 


dx 


(6) 


The  above  system  of  ordinary  differential  equations  has  been  solved  using  the  CHEMKIN 
software  (Kee  et  al.,  1994).  A  set  of  reactions  was  considered  which  included  only  one  electron- 
molecule  collision  reaction  per  each  major  constituent  (N2,  O2,  H2O).  The  dissociation  reactions 
considered  which  were  dependent  upon  the  electric  field  are: 


+e--N-  +N-  +e-  R  =  -k^  [e]  [^^2] 


(7) 
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0^+e--0-+0CD)+e-  R  =  -k^[e][0^ 


(8) 


Hp+e  --^H-  +-OH+e '  R  =  -k^[e\  [Hp] 


(9) 


Reaction  (7)  represents  the  production  of  nitrogen  free  radicals  by  all  possible  mechanisms 
including  dissociation,  ionization,  and  excitation.  The  kinetic  rate  constant  for  this  global  reaction 
is  dependent  upon  the  electric  field  applied  and  hence  would  vary  with  the  external  applied  field. 
Similarly,  the  rate  constants  for  reactions  (8)  and  (9)  represent  the  global  kinetic  rate  constants  for 
the  production  of  oxygen  and  hydroxyl  radicals,  respectively.  Although  ionization  is  a  major  process 
taking  place  in  the  corona  discharge,  it  was  assumed  in  the  present  model  that  the  ions  do  not  play 
a  major  role  in  the  chemistry  of  the  neutral  gas  molecules. 

The  set  of  reactions  used  for  modeling  NO  removal  in  the  case  of  nitrogen,  dry  air,  and  humid 
air  is  shown  in  Table  1 .  The  rate  constants  for  these  reactions  are  obtained  primarily  from  the 
atmospheric  chemistry  literature  (Willis  &  Boyd;  1976,  Baulch  et  al.;  1983,  Demore  et  al.,  1987;  and 
Atkinson  et  al,  1989),  and  were  cross-checked  with  the  NIST  Database  (NIST  Standard  Reference 
Data  Base  17.0,  1994).  The  criterion  for  choosing  the  reactions  was  to  consider  reactions  of  neutral 
species  having  significant  kinetic  rate  constants  (>  1.0  x  10’  cm^molc'sec'). 
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Table  1  ;  List  of  Chemical  Reactions  Considered  for  Modeling 


Chemical  Reactions  Rate  Constants  Source 

[cm^  .mole~^  .sec~^) 


N2  — >  27V- 
O2  — *•  O'  +  0(^  D) 


+  7V-  — 

^N2 

+  0-  — 

^  NO 

+  NO  - 

A^2  +  0- 

+  O2  - 

^  NO  +  O' 

+  NO2  - 

—^2NO 

+  NO2  " 

N2O  +  0 

O'  +  O'  — >  O2 
O'  +  O2  — ^  O3 


IU[E] 

K2[E] 

6.38  X  10^° 
1.45  X  10" 
2.04  X  10^^ 
5.36  X  10' 
1.81  X  10" 
1.81  X  10" 
1.02  X  10" 
7.83  X  10® 


O'  +  O3  — ^  2O2  4.82  X  10® 

0{^D)  — ^  O'  6.95  X  10^ 


0{^D)  +  N02 

— >  NO +  O2 

2.95  X  10" 

0(^D)  +  O3  — 

2O2 

7.23  X  10" 

OCD)  +  O3  - 

->  20-  -h  O2 

7.23  X  10" 

OCD)  +  N2O 

— >  7V2  +  O2 

2.95  X  10" 

0{^D)  +  N2O 

— ^  2NO 

4.04  X  10" 

NO  +  O3  — > 

NO2  +  O2 

8.43  X  10® 

NO  +  0-  — ^ 

N'+02 

2.77  X  10® 

Expt.  Data. 

Expt.  Data. 

Willis  h  Boyd,  1976. 
Willis  k.  Boyd,  1976. 
Willis  k  Boyd,  1976. 
Atkinson  et  al.,  1989. 
.4tkinson  et  al..  1989. 
Atkinson  et  al.,  1989. 
Willis  k  Boyd,  1976. 
Willis  k  Boyd.  1976. 
Willis  k  Boyd.  1976. 
Demore  et  al..  1987. 
Schofield,  1979. 
Atkinson  et  al..  1989. 
.Atkinson  et  al..  1989. 
Demore  et  al.,  1987. 
Demore  et  al.,  1987. 
Demore  et  al.,  1987. 
Demore  et  al.,  1987. 


continued  on  next  page 
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continued  from  previous  page 


Chemical  Reactions  Rate  Constants 

{cm^  ,mole~^  ,sec~^) 

NO  +  O- — >N02  1.40  x  1012 

NO  +  NO3  — ^  2NO2  1.75  X  101^ 

NO  +  NO3 — y2N0  +  02  1.84  x  1011 
N  O2  "f"  O3  — ^  NO3  “I"  O2  1.93  X  10' 

NO2  +  0-  — yNO  +  02  5.84  x  IO12 

NO2  +  NO3  — ^ 

NO  "j"  A  O2  "I"  O2  6.4.5  x  10^ 

.V2O4  +  M  — ^  2NO2  +  M  2.06  X  10^ 

2NO3  2NO2  +  O2  1.92  X  10» 

N03  +  0- — >N02  +  02  1.02  x  1013 

O3  +  N-  — ^  +  O2  6.02  X  10' 

2N O2  "I"  M  — >  N2O4  -j-  AI  .3.06  X  lOi^i 
N02  +  0-^N03  1.29  x  1012 

H2O  —^H-  +  OH-  K3[E] 

N-  +  OH- — >NO  +  H-  1.45x1013 

0(1/))  +  H2—^H-  +  oh-  6.02  X  1013 
0{^D)  +  H20 — ^20H-  1.33  X  lOi'i 

0-+H--^OH-  1.08x1012 

0-  +  HO-2  — ^  OH-  +  O2  3.4  x  1013 
0-  +  OH- — ^02  +  H-  1.99x1013 

OH-  +  H2-^  H2O  +  H-  4.03  X  10^ 

OH-  +  OH-  H2O  +  0-  1.08  X  1012 


Source 


Atkinson  et  al.,  1989. 
Demore  et  al.,  1987. 
Sutherland  et  al.,  1975 
Atkinson  et  al.,  1989. 
Atkinson  et  al.,  1989. 

Baulch  et  al.,  1983. 
Baulch  et  al.,  1983. 
Baulch  et  al.,  1983. 
Demore  et  ah,  1987. 
Baulch  et  al.,  1983. 
Baulch  et  al.,  1987. 
Atkinson  et  ah,  1989. 
Expt.  Data 
Demore  et  ah,  1987. 
Demore  et  ah,  1987. 
Demore  et  al.,  1987. 
Willis  &  Boyd,  1976. 
Atkinson  et  ah,  1989. 
Atkinson  et  ah,  1989. 
Atkinson  et  ah,  1989. 
Atkinson  et  ah,  1989. 
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continued  from  previous  page 
Chemical  Reactions 


OH-  +  OH-  H2O2 

OH-  +  HO-^  H2O  +  O2 

OH-  +  H2O2  H2O  +  HO-^ 

OH-  +  H-  ^  H2O 
H-  +  O2  HO-2 

h-  +  03^  oh-  +  O2 

H-  +  HO,  — >  H2  +  O2 
H-  +  HO:,  — >  -lOH- 
H-  +  H0-,  — ^  H2O  +  0- 
H-  +  H2O2  H2O  +  OH- 
H0-,  +  HO:,  ^  H2O2  +  O2 
NO  +  OH-  HNO2 
NO2  +  OH-  — ^  HNO3 
NO2  +  OH-  — ^  HO-^  +  NO 
HNO2  +  0H-  ^  N02^  H2O 
HNO3  +  OH-  — >  H2O  +  NO3 
OH-  ^  H-  ^  H2  +  O- 
H2O2  ^H-  — >  HO-^  +  H2 
H2O2  +  0-  ^  OH-  +  HO-^ 

O3  +  HO:,  — ^  OH-  +  2O2 
O3  +  OH-  HO-2  +  O2 
HNO2  +  0-  NO2  +  OH- 
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1.81  X  10^ 


continued  from  previous  page 
Chemical  Reactions 


Rate  Constants  Source 
{cm^  .mole~^  .sec~^) 


HNOz  +  HNO2  — ^ 

NO2  +  H2O  9.64  X  10®  Kaiser  et  ah.  1977. 

HNOz  +  0-  ^  NOz  +  OH-  1.81  x  10®  Kaiser  et  al.  1978. 

HNOz  +  H-  — >  HNO2  +  OH-  6.02  x  10^°  Baulch  et  ah.  1983. 

NOz  +  H-  — NO2  +  OH-  6.63  x  10^®  Boodaghians  et  al.,  1988. 

NOz  +  OH-  — »•  NO2  +  H 0-2  1.39  x  10^®  Boodaghians  et  ah,  1988. 

NOz  +  HO-2  — ^  HNOz  +  O2  7.63  x  lO^^  Hoffmann  et  ah,  1988. 

NOz  +  HO2  — ^ 

NO2  +  OH-  +  O2  3.05  X  10^^  Hoffmann  et  ah.  1988. 

N-  +  HO-2  — ^  NO  +  OH-  1.33  x  10^®  Brune  et  al..  1983. 

NO  +  HO-2  NO2  +  OH  4.66  x  lO^^  NIST  Database. 

NO  +  HO-2  — >  HNOz  8.43  x  lO^®  NIST  Database. 
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SECTION  IV 

EXPERIMENTAL  RESULTS  &  DISCUSSION 

A.  PHYSICAL  DISCHARGE  CHARACTERISTICS 

The  physical  aspects  of  the  corona  discharge  including  the  corona  onset  voltage,  the  peak 
voltage,  the  pulse  rise  time,  and  the  pulse  width  were  measured  at  the  location  where  the  pulsed 
discharge  was  delivered  to  the  centrally  located  high  voltage  discharge  wire  electrode. 

First,  a  power  supply  calibration  was  made  to  find  the  relation  between  the  applied  AC 
voltage  delivered  to  the  pulse  forming  network  and  the  voltage  indicated  on  the  control  section  of 
the  transformer/rectifier  (dial  voltage).  A  plot  of  the  applied  voltage  versus  the  dial  voltage  is  shown 
in  Figure  6.  The  applied  voltage  refers  to  the  peak  AC  output  voltage  of  the  transformer  as  measured 
by  the  high  voltage  probe  and  the  oscilloscope.  This  figure  shows  an  apparent  linear  relationship 
between  the  dial  voltage  and  the  applied  voltage.  Linear  regression  was  performed  on  the  data,  and 
the  correlation  coefficient  was  found  to  be  0.9998  and  the  slope  of  the  line  to  be  0.964. 


Figure  6.  Power  Supply  Calibration. 

Pulsed  corona  onset  (corona  initiation  voltage)  in  the  reactor  was  found  to  occur  at  a  pulsed 
input  voltage  of  between  13  and  14  kV  in  an  atmosphere  of  dry  air  and  at  about  10  kV  in  a  dry 
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nitrogen  atmosphere.  Below  the  onset  of  corona,  the  voltage  and  the  current  readings  measured  at 
the  reactor  using  the  oscilloscope  were  very  small.  Above  the  corona  onset  the  peak  voltage  and  the 
current  increased  dramatically.  The  rise  time  of  the  voltage  pulse  was  found  to  vary  between  1 5  and 
30  nanoseconds. 

A  plot  of  a  typical  voltage  waveform  is  shown  in  Figure  7.  As  seen  in  the  figure,  there  is  a 
sharp  increase  in  the  voltage,  i.e.,  a  very  short  rise  time,  followed  by  an  exponential  decay  of  the 
voltage.  The  pulse  width  at  half-maximum  of  the  voltage  pulse  varied  from  500  to  700  nanoseconds. 
The  current  waveform  is  shown  in  Figure  8.  The  current  pulses  were  shorter  in  width  (<  300 
nanoseconds)  when  compared  to  the  voltage  pulses.  The  power  input  per  pulse  was  calculated  by 
multiplying  the  voltage  and  the  current  waveforms  and  measuring  the  area  under  the  resulting 
waveform.  A  plot  of  the  pulsed  power  waveform  is  shown  in  Figure  9.  Table  2  summarizes  the 
peak  voltage,  the  rise  time,  and  the  pulse  width  of  the  pulsed  voltage  as  a  function  of  dial  voltage  in 
a  dry  air  atmosphere. 

As  the  table  shows,  the  peak  pulsed  voltage  was  significantly  higher  than  the  dial  voltage, 
and  both  the  rise  time  and  the  pulse  width  increased  with  an  increase  in  the  applied  voltage. 

As  originally  conceived,  the  pulse  generating  circuit  did  not  include  a  load  resistor  in  parallel 
with  the  reactor  electrodes  (see  Figure  3).  However,  it  was  found  that  the  spark  gap  would  not 
fire  even  at  high  input  voltages  because  the  effective  resistance  of  the  gap  between  the  reactor 


Table  2:  Voltage  Characteristics  of  the  Pulse  in  Dry  Air 


Dial  Voltage 

Peak  Voltage 

Rise  Time 

Pulse  Width 

(kV) 

(kV) 

(ns) 

(ns) 

30 

38.8 

13.5 

449 

35 

44.8 

16.5 

616 

40 

51.4 

18.0 

627 

50 

59.6 

29.0 

694 
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Figure  9  Power  Waveform  in  an  Atmosphere  of  Dry  Air  at  40  kV 
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electrodes  was  too  high.  In  order  to  decrease  the  effective  resistance  of  the  gap  between  the  outer 
grounded  cylinder  and  the  high  voltage  wire  electrode,  a  load  resistance  was  connected  from  the 
spark  gap  output  to  ground  in  parallel  to  the  two  reactor  electrodes.  The  nature  of  the  voltage 
waveform  was  affected  by  the  value  of  the  load  resistance.  As  the  resistance  of  the  load  resistor  was 
decreased  from  300  Q  to  77.1  Q,  the  decay  time  (time  required  for  the  voltage  to  drop  from  its 
maximum  value  to  its  minimum  value)  for  the  voltage  pulses  increased.  A  load  resistor  of  300  Q 
was  used  in  all  the  experiments  conducted  for  NO  removal. 

In  this  study,  the  pulse  frequency  was  kept  constant  at  60  Hz  and  the  voltage  was  varied  in 
order  to  increase  the  energy  input  into  the  reactor.  Table  3  shows  the  peak  value  of  the  current 
waveform  and  the  corresponding  energy  per  unit  pulse  calculated  as  a  function  of  applied  voltage. 
The  power  delivered  to  the  reactor  was  determined  by  multiplying  the  energy  per  pulse  by  the 
frequency  of  pulsing.  With  the  Tektronix  current  probe  P6021,  the  energy  input  varied  from  8.0 
mJ/pulse  to  148.0  mJ/pulse  in  the  case  of  dry  air.  No  differences  in  the  pulse  characteristics  were 
observed  when  humid  air  was  fed  to  the  reactor. 


Table  3:  Current  &  Power  Characteristics  of  the  Pulse  in  Dry  Air 


Dial  Voltage 

(kV) 

Peak  Current 

(A) 

Energy  per  Pulse 

(mJ) 

Power 

(watts) 

30 

42 

8.0 

0.48 

35 

52 

20. 

1.20 

40 

64 

56. 

3.36  - 

50 

82 

148. 

8.88 

The  energy  delivered  per  pulse  was  close  to  that  reported  by  Tas  (1995)  (6-7  mJ/pulse).  The 
power  was  found  to  vary  from  0.48  to  8.88  watts.  Mizuno  et  al.  (1995)  reported  an  input  power  of 
20  to  40  watts  which  is  higher  than  the  values  obtained  in  the  present  study.  The  power  density 
varied  from  1 .4  to  158  J/1  which  is  in  the  range  of  that  reported  by  Wolf  et  al.  (1996)  and  Penetrante 
etal.  (1996). 

Experiments  were  designed  to  determine  the  extent  of  NO  removal  from  several  gases 
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with  different  compositions  in  order  to  study  the  reaction  mechanisms  for  NO  removal.  The 
experiments  were  conducted  at  four  different  residence  times  and  at  four  different  voltages.  All 
of  the  experiments  were  repeated  twice  to  insure  the  reproducibility  of  the  results.  The  average 
concentration  determined  in  two  similar  runs  was  used  to  plot  the  figures  of  NO  and  NO2 
concentrations  as  functions  of  residence  time  and  operating  voltages.  As  an  initial  case, 
experiments  on  NO  removal  in  an  atmosphere  of  nitrogen  was  considered.  This  was  followed  by 
experiments  in  dry  air,  humid  air,  and  finally  in  the  presence  of  gas  additives  including  carbon 
monoxide  and  ethylene. 

B.  NO  REMOVAL  IN  NITROGEN 

NO  removal  experiments  in  an  atmosphere  of  dry  nitrogen  were  conducted  at  different 
residence  times  (7.2, 14.5,  29.0,  and  44.0  seconds)  and  at  different  dial  voltages  (30,  35, 40,  and 
50  kV).  The  initial  concentration  of  NO  during  the  experiments  was  100  ppmv. 

The  NO  removal  data  as  a  function  of  residence  time  for  different  voltages  is  shown  in 
Figure  10a.  It  was  observed  that  the  NO  concentration  decreased  as  the  residence  time  in  the 
reactor  was  increased  at  a  particular  voltage.  In  addition,  as  the  applied  voltage  was  increased  for 
a  particular  residence  time,  the  NO  concentration  decreased.  Complete  removal  of  NO  was 
observed  at  the  higher  voltages  of  40  and  50  kV  at  a  residence  time  of  44.0  seconds. 

The  NO,  concentration  profile  as  a  function  of  residence  time  for  different  dial  voltages  is 
shown  in  Figure  10b.  A  maximum  of  8  ppm  NO2  formation  was  measured  with  the  NO,,  monitor 
at  the  voltages  of  40  and  50  kV.  The  concentration  profiles  of  NO2  at  35,  40,  and  50  kV  passed 
through  maxima  as  the  residence  times  in  the  reactor  increased.  Also,  as  the  voltage  was 
increased,  the  time  corresponding  to  the  maximum  concentration  decreased.  At  30  kV  the  NO, 
concentration  was  steadily  increasing  with  the  residence  time.  These  profiles  indicate  that  at 
higher  voltages  (35,  40.  &  50  kV),  NO2  is  an  intermediate  product,  and  at  longer  residence  times 
NO,  was  reduced  back  to  molecular  nitrogen  and  oxygen.  At  30  kV,  the  amount  of  power  input 
was  not  sufficient  to  reduce  the  NO,  formed.  Approximately  97%  removal  of  NO,  was  obtained 
at  high  voltages  (40  and  50  kV)  and  at  a  residence  time  of  44.0  seconds. 

Penetrante  et  al.  (1996)  and  Wolf  et  al.  (1996)  plot  the  NO  concentration  versus  the 
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power  density.  Figure  1 1  shows  our  NO  and  NO2  data  ploted  versus  power  density  with  a 
comparison  to  the  theoretical  predictions  of  Penetrante.  The  data  of  Penetrante  et  al.  (1996)  and 
Wolf  et  al.  (1996)  for  lower  voltages  and  applied  fields  follow  the  theoretical  line  of  Penetrante 
et  al.  (1996),  however,  the  data  of  Wolf  et  al.  (1996)  for  higher  field  intensities  and  our  data 
show  a  significantly  lower  power  density  for  a  given  NO  removal. 

The  kinetic  model  described  earlier  was  used  to  determine  the  rate  constant  for  the 
nitrogen  dissociation  reaction  by  fitting  the  model  to  the  NO  experimental  data.  The  number  of 
chemical  reactions  considered  to  calculate  the  concentration  profiles  of  different  species  was  31. 
Figure  12  shows  the  mole  fi-actions  of  NO,  NO2,  and  N2O  calculated  by  the  model  as  a  function 
of  residence  time  at  50  kV.  Experimental  values  of  the  NO2  concentration  profile  at  50  kV  in 
Figure  1 1  qualitatively  match  with  the  model  profile  of  NO2  at  50  kV  shown  in  Figure  12  (100 
ppm  =  0.0001  in  terms  of  mole  fraction).  The  maximum  value  of  NO2  concentration  in  the 
model  profile  is  around  17  ppm  whereas  the  maximum  value  found  by  the  experiments  at  50  kV 
was  around  8  ppm.  Figure  12  also  shows  a  small  formation  of  N2O.  Essentially  the  model 
predicts  around  90%  reduction  of  initial  NO  to  molecular  nitrogen  and  oxygen  in  an  atmosphere 
of  nitrogen  at  50  kV. 

The  model  profile  of  NO  and  the  experimental  data  are  shown  in  Figure  13.  As  seen  in 
Figure  13,  the  model  profiles  match  very  well  with  the  experimental  values  at  35,  40,  and  50  kV 
at  longer  residence  time.  There  is  a  slight  deviation  of  the  model  from  the  experimental  points  at 
30  kV.  The  experimental  data  of  NO  removal  at  7.2  seconds  residence  time  for  all  the  voltages 
except  50  kV  were  less  than  that  predicted  by  the  model. 

The  rate  constant  of  nitrogen  dissociation  as  a  function  of  applied  voltage  thus  obtained 
by  fitting  the  model  to  the  experimental  data  is  shown  in  Figure  14.  This  figure  shows  an  almost 
linear  relation  of  the  rate  constant  as  a  function  of  applied  electrical  field.  Of  course,  since  the 
field  is  non-linear  with  power  as  shown  in  Table  3,  the  reaction  rate  constant  is  a  non-linear 
function  of  the  power  delivered  to  the  reactor.  Many  theoretical  studies  of  pulsed  corona  and 
other  non-thermal  plasma  reactors  report  the  pseudo  first-order  rate  constant  for  electron  collison 
reaction  to  be  of  the  form  k -[E/n]k,  where  k'  is  the  pseudo  first  order  rate  constant,  [E/n]  is  the 
reduced  field,  i.e.  value  of  electrical  field  divided  by  the  gas  density,  and  k  is  the  true  second- 
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Mole  Fraclion  o(  NO  Wo'®  Fraclion 


Figure  12  Model  Profile  of  Different  Species  in  an  Atmosphere  of  Nitrogen 


Figure  13  Model  Fit  for  NO  Concentration  Profile  in  an  Atmosphere  of  Nitrogen 


order  rate  constant  for  the  electron  molecule  collision.  For  constant  gas  density,  this  relationship 
would  predict  a  linear  dependence  upon  the  electrical  field  as  seen  in  our  study. 

The  rate  constant  determined  is  a  pseudo  first-order  rate  constant,  which  is  essentially  the 
electron  concentration  multiplied  by  the  actual  second-order  dissociation  rate  constant  as  in 
Equations  7,  8,  and  9.  Both  experimental  and  modeling  results  indicate  that  the  main  reaction 
pathway  for  NO  removal  in  nitrogen  is  through  reduction  to  molecular  nitrogen.  The  main 
reactions  responsible  for  the  removal  are: 


N^+e  ^  N-+N-+e 


(10) 


N0*0-  -  NO^ 


(11) 


NO+N-  -  N^+0- 


(12) 


NO^^N-  -  Np+0-  (13) 

NO^+N-  -  A/’2+6>2  (14) 

As  the  applied  voltage  is  increased,  the  amount  of  nitrogen  radicals  produced  by  Equation 
10  increases.  The  reactions  responsible  for  NO  removal  are  Equations  1 1  and  12.  Equation  1 1 
results  in  the  formation  of  NO,  but  Equations  13  and  14  are  fast  enough  to  reduce  the  NO, 
formed  back  to  N2  and  to  N,0. 

Tokunaga  et  al.  (1984)  used  an  electron  beam  to  consider  NO  removal  as  a  function 
of  dose  rate  in  an  atmosphere  of  nitrogen,  and  concluded  that  NO  was  mainly  decomposed  to 
nitrogen  and  oxygen.  The  products  from  the  irradiation  were  analyzed  and  NO,  and  H  0  were 
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Figure  14.  Rate  Constants  for  Nitrogen  Dissociation  in  Pure  Nitrogen 

present  in  very  small  amounts  (NO2=10%,  N2  0=3%).  Tas  (1995)  used  pulsed  streamer  corona 
discharge  to  consider  NO  removal  in  an  atmosphere  of  nitrogen,  and  concluded  that  the  reduction 
pathway  was  the  main  path  for  the  removal  of  NO  in  nitrogen.  This  conclusion  is  also  supported 
by  other  studies  (Wolf  et  al.,  1996;  Penetrante  et  ah,  1996). 

Penetrante  et  al.  (1995)  considered  in  detail  the  kinetics  of  corona  discharge  by  solving  for 
the  electron  energy  distribution  function  to  calculate  the  kinetic  rate  constants  for  the  electron-gas 
reactions  in  order  to  predict  the  profile  of  different  species  for  the  case  of  NO  removal  in  an 
atmosphere  of  nitrogen.  The  chemical  kinetics  calculation  included  287  reactions  and  37  species 
involving  charged  as  well  as  neutral  species.  This  model  also  predicts  the  reduction  of  NO  as  the 
main  pathway  for  NO  removal.  It  is  important,  however,  to  note  that  in  the  present  study  the  model 
developed  for  the  case  of  NO  removal  in  nitrogen  has  considered  only  13  neutral  species  and  31 
chemical  reactions  as  compared  to  the  vast  number  of  reactions  considered  by  Pentrante  et  al.  (1995). 
It  is,  therefore,  possible  to  characterize  the  chemical  reactions  taking  place  in  the  pulsed  streamer 
corona  reactor  by  considering  the  neutral  species  reactions. 
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C.  NO  REMOVAL  IN  DRY  AIR 

Experiments  on  NO  removal  in  an  atmosphere  of  dry  air  at  various  voltages  (30,  35, 40,  and 
50  kV)  and  at  different  gas  residence  times  (3.6, 7.2, 14.5, 29.0,  and  60.0)  were  conducted  to  study 
the  effect  of  oxygen  addition  on  the  NO  and  NO2  concentration  profiles. 

The  feed  air  used  for  the  study  was  from  a  compressed  air  source  from  the  building.  The 
compressed  air  source  was  passed  through  a  series  of  drier  columns  (anhydrous  calcium  sulfate)  to 
remove  the  water  vapor  in  the  air.  The  initial  NO  concentration  was  about  116  ppm. 

The  concentration  of  NO2  varied  from  0  ppm  to  25  ppm  as  the  residence  time  in  the  reactor 
(without  corona)  increased.  This  indicated  that  a  small  fraction  of  the  initial  NO  was  oxidized  to 
NO2  when  mixed  with  air  in  the  reactor  even  in  the  absence  of  corona  discharge.  The  temperature 
of  the  gas  varied  from  12.9  to  19  "C  during  the  experiments. 

NO  concentration  profiles  as  a  function  of  residence  time  for  different  pulsed  voltages  are 
shown  in  Figure  15.  NO  concentration  decreased  as  the  applied  voltage  was  increased  for  a 
particular  residence  time.  It  also  decreased  as  the  residence  time  was  increased  for  a  particular 
voltage.  Comparing  Figures  10  and  15,  the  rate  of  NO  removal  was  faster  in  the  case  of  dry  air  than 
in  nitrogen.  At  50  kV  in  an  atmosphere  of  nitrogen  90%  of  NO  was  removed  whereas  in  dry  air 
almost  100%  removal  was  observed  at  14.5  seconds  residence  time.  Using  an  initial  concentration 
of  250  ppmv  of  NO,  Chang  et  al.  (1992)  observed  that  at  room  temperature,  the  addition  of  oxygen 
reduces  the  NO  removal  rate  in  case  of  dielectric  barrier  discharge.  This  is  in  contradiction  to  what 
was  observed  in  the  experiments  conducted  in  this  study.  Tas  (1995)  also  observed  that  for  an  initial 
concentration  of  1200  ppm,  the  NO  removal  decreased  with  an  increase  in  the  oxygen  composition 
on  the  feed  gas. 

The  NOj  concentration  profiles  are  shown  in  Figure  16.  The  NO2  concentration  increased 
with  residence  time  for  all  voltages  at  short  residence  times.  The  NO;  profiles  at  40  and  50  kV  go 
through  maxima  and  the  concentrations  decrease  slightly  at  longer  residence  times.  NO^ 
concentration  profiles  at  30  and  35  kV  increase  steadily  with  residence  time.  Unlike  the  NO2  profile 
in  the  case  of  NO  removal  in  nitrogen,  the  removal  of  NO2  in  dry  air  was  not  observed  at  longer 
residence  times  and  at  high  voltages. 

The  results  indicate  that  most  of  the  NO  was  converted  into  NO2  at  the  tested  voltages  and 
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Figure  15  NO  Concentration  Profile  in  an  Atmosphere  of  Dry  Air 


Figure  16  NO2  Concentration  Profile  in  an  Atmosphere  of  Dry  Air 


residence  times.  Total  NO^  was  higher  than  the  initial  NO  concentration  at  14.5  and  29.0  seconds 
residence  time  and  at  50  kV.  This  indicated  that  there  was  NO^  formation  from  the  dry  air  itself 
during  the  discharge.  To  confirm  this,  an  experiment  to  consider  the  formation  of  NO  and  NOj  in 
the  pulsed  corona  discharge  in  dry  air  was  conducted.  It  was  observed  that  at  60  seconds  residence 
time  and  at  an  applied  voltage  of  50  kV,  the  formation  ofNOj  was  10  ppm.  Tokunaga  et  al.  (1984) 
also  noticed  an  increase  in  the  total  NO,;  concentration  when  the  initial  concentration  of  NO  was 
around  50-100  ppm.  For  an  initial  concentration  of  100  ppm  of  NO,  the  NO^  concentration  at  a  dose 
rate  of  5  Mrad  was  120  ppm,  and  it  was  found  to  increase  with  dose  rate  for  a  gas  composition  of 
3%  oxygen  in  an  atmosphere  of  nitrogen. 

A  measurement  of  N2O  was  made  for  the  44.0  second  residence  time  experiment  at  50  kV, 
and  the  concentration  was  found  to  be  17  ppm.  The  measurement  was  made  by  a  GC/MS  analysis 
of  the  gas  sample  coming  out  of  the  outlet  of  the  reactor. 

For  the  case  of  NO  removal  in  dry  air,  the  two  unknown  rate  constants  corresponding  to  the 
nitrogen  and  oxygen  dissociation  were  obtained  by  fitting  the  model  profiles  with  the  experimental 
data  for  NO  and  NO,.  The  model  profiles  of  NO,  NOj,  N,0,  and  0-  are  shown  in  Figure  1 7  for  the 
treatment  cases  at  50  kV.  NO  removal  corresponded  with  NO,  formation  at  shorter  residence  times, 
and  at  the  longer  residence  times  the  NO,  concentration  leveled  out  at  80%  of  the  initial  NO.  The 
O3  concentration  increased  only  after  a  significant  lag  time  where  the  NO  concentration  decreased. 
This  is  due  to  the  fact  that  in  the  presence  of  NO  the  amount  of  ozone  formation  due  to  oxygen 
radicals  reacting  with  oxygen  molecules  is  very  small  when  compared  to  oxygen  radicals  reacting 
with  NO  to  produce  NO2.  Once  the  NO  molecules  were  removed,  ozone  formation  was  unhindered 
and  the  slope  of  the  ozone  concentration  as  a  fimction  of  residence  time  started  increasing.  The  NjO 
concentration  as  predicted  by  the  model  for  a  residence  time  of  44.0  seconds  at  50  kV  is  around  15 
ppm  which  matches  well  with  the  experimentally  observed  value  .  This  is  a  significant  observation 
since  N,0  data  was  not  used  in  fitting  the  model,  i.e.,  the  N  p  measurement  provides  an  independent 
check  of  the  model.  In  case  of  NO,,  however,  there  was  around  20%  deviation  between  the  model 
profile  and  the  experimental  values  at  50  kV. 

The  model  profile  for  NO  and  the  experimental  points  used  to  fit  the  model  are  shown  in 
Figure  1 8.  As  seen  in  the  figure,  there  is  a  deviation  between  the  experimental  points  and  the  model 
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Figure  17  Model  Profile  of  DiiTerent  Species  in  an  Atmosphere  of  Dry  Air 


Figure  18  Model  Fit  for  NO  Concentration  Profile  in  an  Atmosphere  of  Dry  Air 


at  short  residence  time  (3.6  seconds).  One  reason  for  this  behavior  could  be  that  at  short  residence 
time,  the  flow  pattern  in  the  reactor  was  highly  turbulent  and  hence  the  reactor  would  operate  more 
like  a  well-mixed  reactor  than  as  a  plug-flow  reactor.  The  model  profile  for  NO  fits  well  with  the 
experimental  points  at  longer  residence  times. 

The  main  reaction  mechanisms  for  the  removal  of  NO  in  an  atmosphere  of  dry  air  was 
through  the  oxidation  of  NO  to  NOj  and  the  main  reactions  in  addition  to  Reactions  1 1  to  14  are: 


Oj+e'  -  0-+0(‘£>)+e-  (15) 

-  O3  (16) 

Reaction  15  produces  the  necessary  oxygen  radicals  required  for  the  oxidation  of  NO. 
Reaction  1 1  converts  the  NO  to  NO,  and  Reaction  13  converts  the  NO  2  formed  into  N  0.  The  ozone 
formation  was  significant  once  the  NO  molecules  are  removed  and  Reaction  1 6  is  the  corresponding 
reaction  pathway. 

The  rate  constants  obtained  for  nitrogen  and  oxygen  dissociation  as  a  function  of  applied  voltage 
are  shown  in  Figure  19.  The  slope  for  the  oxygen  dissociation  is  higher  than  that  for  nitrogen 
dissociation  in  air.  The  nitrogen  dissociation  rate  constant  in  dry  air  was  not  the  same  as  in  the  case 
of  a  pure  nitrogen  feed  gas.  The  rate  constant  was  0.0475  s'  in  an  atmosphere  of  nitrogen  at  50  kV 
but  in  dry  air  the  rate  constant  was  reduced  to  0.03  s’'  at  the  corresponding  voltage.  The  presence 
of  electronegative  gas  such  as  oxygen  is  known  to  have  an  effect  on  the  nitrogen  dissociation  (Chang 
et  al.,  1992). 

The  differences  in  the  rate  constants  for  nitrogen  and  oxygen  dissociation  in  dry  air  are 
functions  of  the  applied  voltage.  At  low  voltage  (30  kV)  the  rate  constants  are  almost  equal 
(approximately  .01  s')  and  at  the  high  voltage  (50  kV)  the  rate  constants  are  0.025  s’'  for  nitrogen 
dissociation  and  0.25  s''  for  oxygen  dissociation.  This  factor  of  ten  difference  is  significantly  smaller 
tlian  theoretical  studies  (Penetrante  et  al.,  1995)  where  it  is  predicted  tliat  tlie  G-values  for  oxygen 
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Figure  19  Rate  Constants  of  Nitrogen  and  Oxygen  Dissociation  as  a  Function  of 
Applied  Voltage  in  an  Atmosphere  of  Dry  Air 


Figure  20  O3  Concentration  Profile  Predicted  by  the  Model 


and  nitrogen  dissociation  would  differ  by  about  600  times.  It  should  also  be  noted  that  these 
theoretical  studies  also  do  not  predict  a  field  dependence  as  seen  in  our  work. 

Tokunaga  et  al.  (1984)  conducted  experiments  for  NO  removal  in  a  dry  mixture  of  nitrogen 
and  oxygen  (N2  97%,  O2  3%)  using  an  electron  beam,  and  they  concluded  that  for  high  initial 
concentrations  of  NO  and  at  high  dose  rate,  the  ratio  of  the  final  concentrations  of  NO  and  NO2  after 
the  treatment  were  given  by  the  empirical  relation  of  Equation  1.  Tas  (1995)  also  considered  the 
removal  of  NO  from  dry  air  and  concluded  that  the  empirical  relation  holds  true.  However, 
experiments  and  modeling  conducted  in  the  present  work  do  not  follow  the  empirical  relation 
suggested  by  Tokunaga  et  al  (1984).  That  relation  holds  true  only  for  cases  in  which  the  initial 
concentration  of  NO  was  very  high  (500-1000  ppm).  Penetrante  (1993)  considered  modeling  the 
concentration  profiles  of  NO  and  NO2  in  an  atmosphere  of  dry  air  at  room  temperature,  and  he 
concluded  that  for  an  initial  concentration  of 400  ppm  of  NO2,  50%  of  the  initial  NO  was  converted 
into  NO2,  which  was  later  again  converted  to  N2O  and  Nj.  In  that  study,  the  model  prediction  for 
an  initial  concentration  of  100  ppm  is  80%  conversion  of  the  initial  NO  to  NO2  and  20  ppm  N2O 
formation. 

D.  OZONE  FORMATION  IN  DRY  AIR 

The  formation  of  ozone  in  dry  air  with  pulsed  streamer  corona  treatment  was  considered  at 
50  kV  and  with  a  residence  time  of  44.0  seconds.  Measurements  by  the  KJ  solution  method  gave 
ozone  concentration  under  these  conditions  of  approximately  670  ppm.  The  amount  of  ozone 
predicted  by  the  model  is  shown  in  Figure  20.  The  model  predicts  200  ppmv  of  ozone  to  be 
produced  at  a  gas  residence  time  of  44.0  seconds  and  at  50  kV,  whereas  the  value  observed  in  the 
experiment  was  670  ppm.  Further  refinement  of  the  model  to  account  for  this  variation  is  under 
development  and  more  extensive  measurements  of  ozone  concentration  under  a  wider  range  of 
conditions  is  being  performed.  It  should  be  noted  that  the  KI  solution  gives  a  fairly  accurate 
measurement  of  ozone  level  and  that  the  thesis  of  Mr.  Kalyana  reported  measurements  of  up  to  1200 
parts  per  million  as  measured  by  a  commercial  ozone  monitor.  This  level  of  ozone  is  too  low  to  be 
accurately  measured  by  the  ozone  monitor  according  to  the  manufacturer's  specifications. 


E.  NO  REMOVAL  IN  HUMID  AIR 

Water  vapor  is  present  in  concentrations  of  up  to  5%  by  mole  in  the  exhaust  gases  produced 
by  the  combustion  of  fossil  fuels  or  from  jet  cell  engines  (Spicer  et  al,  1992).  The  hydroxyl  radicals 
produced  from  the  pulsed  corona-induced  dissociation  of  water  molecules  react  with  NOj  to  form 

HNO3. 

Experiments  were  conducted  in  order  to  test  the  removal  of  NO  in  humid  air.  The  feed  air 
was  humidified  by  passing  it  through  a  humidifier  bringing  the  water  content  to  83%  relative 
humidity  at  18  "C.  This  corresponds  to  1.5%  of  gas  composition  by  volume.  As  in  previous 
experiments  the  voltage  was  varied  at  different  residence  times  (7.2, 14.5, 29.0,  and  44.0  seconds), 
and  the  NO  and  NO,  outlet  concentrations  from  the  reactor  were  monitored. 

Figure  21  shows  the  concentration  profile  of  NO  as  a  function  of  residence  time  for  different 
dial  voltages.  When  compared  with  Figure  15  corresponding  to  NO  removal  in  dry  air,  the  amount 
of  NO  removal  for  the  operating  voltages  of  35, 40,  and  50  kV  in  humid  air  was  less  than  that  of  dry 
air.  At  an  operating  voltage  of  50  kV,  80%  removal  of  NO  was  observed  at  a  residence  time  of  7.2 
seconds  in  dry  air.  In  the  presence  of  water  at  the  same  operating  conditions,  75%  removal  of  NO 
was  observed.  Hence  the  addition  of  water  did  not  enhance  the  NO  removal.  However,  with 
dielectric  discharge,  Chang  et  al.  (1992)  reported  an  increase  in  the  removal  efficiency  of  NO  with 
the  addition  of  water  for  an  initial  concentration  of  250  ppmv  of  NO. 

Figure  22  shows  the  concentration  profiles  of  NOj  as  functions  of  residence  time  and 
operating  dial  voltage.  The  NO,  concentration  undergoes  a  maximum  at  about  14.5  seconds 
residence  time  for  35  and  40  kV.  The  maximum  for  50  kV  was  at  7.2  seconds  while  the  30  kV 
profile  had  a  steady  increase  in  the  concentration.  The  concentration  of  NO2  was  less  than  5  ppm 
for  all  the  voltages  at  a  residence  time  of  44.0  seconds.  The  maximum  concentration  of  NO.  attained 
was  22  ppm  for  tire  case  of  50  kV.  Around  95%  removal  of  total  NO,  was  observed  for  40  and  50 
kV  voltages  at  a  gas  residence  time  of  44.0  seconds.  Thus,  with  the  addition  of  water,  the  NO, 
concentration  decreased  significantly.  The  NOj  concentration  was  2.5  ppm  at  50  kV  at  a  residence 
time  of  29.0  seconds  in  the  case  of  humid  air  compared  to  86.4  ppm  in  the  case  of  dry  air  at 
corresponding  conditions. 

In  order  to  measure  the  possible  formation  of  HNO3  from  NO,,  the  exhaust  gases  from  the 
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NO  Concentration  Profile  in  an  Atmosphere  of  Humid  Air 


Figure  22  ;  NO2  Concentration  Profile  in  an  Atmosphere  of  Humid  Air 


reactor  were  sparged  through  de-ionized  water.  The  operating  voltage  was  50  kV  and  the  gas 
residence  time  was  44.0  seconds.  The  liquid  sample  was  then  analyzed  for  nitrate  by  ion 
chromatography.  The  concentration  of  NOj''  in  the  sample  was  900  ppb  by  weight  which 
corresponds  to  1 .64  ppmv  of  HNO3  in  the  gas  sample.  The  amount  measured  as  nitrates  is 
significantly  lower  when  compared  with  the  removal  of  NO2. 

The  important  reactions  which  are  responsible  for  the  reduced  concentration  of  NO,  are: 


H^O+e'  -  H-+OH-+e- 

(17) 

NO+OH-  -  HNO^ 

(18) 

NO^*OH-  -  HNO^ 

(19) 

The  model  profiles  of  NO,  NOj,  HNO2,  and  HNO3  are  shown  in  Figure  23  corresponding  to 
pulsed  voltages  of  50  kV.  The  number  of  chemical  reactions  considered  to  model  the  reactions 
taking  place  in  the  presence  of  water  in  the  system  was  74.  The  NO2  profile  undergoes  a  maximum 
concentration  of  50  ppm  at  16  seconds  residence  time.  The  model  predicts  up  to  35%  conversion 
of  initial  NO  to  HNO3.  HNO  2  is  an  important  intermediate  product  and  is  an  important  route  for  NO 
removal.  The  nitric  acid  produced  tends  to  form  aerosols.  Future  work  should  consider  the  study 
of  aerosol  formation. 

Figure  24  shows  the  NO  model  profile  versus  the  experimental  data  used  to  fit  the  model  for 
calculating  the  rate  constant  for  water  dissociation.  As  the  graph  shows,  the  model  matches  the 
experimental  points  at  short  residence  times,  but  deviates  from  it  at  longer  residence  times.  In 
calculating  the  rate  constant  for  the  water  molecule  dissociation,  an  assumption  has  been  made  that 
the  rate  constants  for  oxygen  and  nitrogen  dissociation  were  not  affected  by  the  addition  of  water. 
That  may  possibly  explain  some  of  the  deviation  of  the  model  from  the  experimental  points. 

The  maximum  concentration  of  NOj  as  observed  in  the  experiments  was  around  20  ppm  for 
the  case  of  50  kV  pulsed  voltage  while  the  model  predicts  a  maximum  of  50  ppm.  The  NO2 
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Figure  23  Model  Profile  of  Different  Species  in  an  Atmosphere  of  Humid  Air 


Figure  24  Model  Fit  for  NO  Removal  in  an  Atmosphere  of  Humid  Air 


concentration  profile  predicted  by  the  model  for  various  voltages  is  shown  in  Figure  25.  From 
Figure  25,  the  NO,  concentration  for  the  cases  of  40  and  50  kV  goes  through  a  maximum,  while  in 
the  cases  of  30  and  35  kV  there  is  a  steady  increase  in  the  concentration  of  NOj.  But  the 
experimental  data  shown  in  Figure  22  indicate  that  at  all  of  the  voltages  except  30  kV,  the  NO^ 
concentration  goes  through  a  maximum. 

The  kinetic  rate  constant  for  the  water  dissociation  obtained  by  fitting  the  model  with  the 
experimental  data  is  shown  in  Figure  26  as  a  function  of  operating  dial  voltage.  The  rate  constant 
increases  monotonically  with  the  applied  voltage.  Tokunaga  et  al.  (1984),  Chang  et  al.  (1992),  and 
Tas  (1995)  observed  that  the  addition  of  water  enhances  the  NO  removal  as  well.  But  in  this  study 
the  addition  of  water  did  not  enhance  the  NO  removal.  One  of  the  reasons  is  that  at  high 
concentrations  of  initial  NO  and  NO,  (~  1000  ppm)  in  dry  air,  the  final  concentrations  of  NO  and 
NO,  come  to  equilibrium  as  suggested  by  the  empirical  relation  given  by  Equation  1.  With  the 
addition  of  water,  NO,  removal  takes  place  due  to  the  formation  of  HNO,  from  NO,.  Since  there 
is  a  decrease  in  the  NO,  level,  the  corresponding  equilibrium  concentration  of  NO  also  decreases 
according  to  the  Equation  1.  Hence  the  enhanced  removal  of  NO  noticed  in  the  presence  of  water 
is  not  due  to  the  direct  reaction  of  NO  with  the  hydroxyl  radicals.  However,  the  experiments 
conducted  in  this  study  indicate  that  the  empirical  Equation  1  does  not  hold  when  the  initial 
concentration  of  NO  is  low  (~100  ppm). 

F.  THE  EFFECT  OF  CO  ON  NO  REMOVAL  IN  HUMID  AIR 

Carbon  monoxide  is  one  of  the  important  constituents  of  combustion  gases.  The  addition 
of  500  ppm  of  CO  to  humid  air  was  studied  to  consider  NO  removal  under  these  conditions. 
Tokunaga  et  al.  (1984)  showed  that  the  addition  of  CO  enhanced  NO  removal  due  to  the  production 
of  HO,  radicals  by  the  following  reaction  mechanism: 


CO^OH-  -  CO^+H- 


(20) 
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H+OH-  -  HO^- 


(21) 


NO+HO-^  -  NO^+OH- 


(22) 


Expcrirncnts  to  dctcrniinc  NO  removal  with  the  addition  of  CO  were  performed  at  different 
voltages  (30,  35,  40,  and  50  kV)  and  at  different  gas  residence  times  (7.2,  14.5,  29.0,  and  44.0 
seconds).  The  inlet  concentration  of  CO  was  maintained  at  500  ppm.  The  concentration  of  NO  and 
NO,  as  a  function  of  residence  time  for  different  operating  voltages  are  shown  in  Figures  27  and  28, 
respectively.  Comparing  the  concentration  profiles  of  NO  removal  without  CO  (Figure  21)  and  with 
CO  (Figure  27),  CO  did  not  have  any  effect  on  the  NO  removal.  The  NO  concentrations  at  7.2 
seconds  residence  time  were  62.6  ppm  and  29.5  ppm  at  40  and  50  kV,  respectively,  in  the  presence 
of  CO.  In  the  absence  of  CO  in  the  system  (see  Figure  21)  the  NO  concentration  at  these  voltages 

for  the  same  residence  time  were  62.0  and  25.6  ppm. 

Figures  22  and  28  can  be  compared  for  NO,  concentration  profiles  to  conclude  that  the 
presence  of  CO  at  500  ppm  does  not  affect  the  NO,  profile.  The  NO,  concentrations  at  7.2  seconds 
residence  time  were  16.6  and  24.9  ppm  for  40  and  50  kV,  respectively,  with  the  presence  of  CO  in 
the  system.  The  corresponding  NO,  concentrations  in  the  absence  of  CO  in  the  system  at  7.2 
seconds  residence  time  was  1 4.8  and  21.8  ppm  for  the  cases  of  40  and  50  kV.  The  amount  of  CO 
added  might  be  too  small  to  provide  a  significant  increase  in  the  NO  removal.  The  minimum 
concentration  with  which  Tokunaga  et  al.  (1984)  observed  any  change  in  the  NO  removal  was  0.29% 
(2900  ppm)  of  CO. 

The  kinetic  model  profiles  of  NO.  NO,,  and  HNO3  for  NO  removal  in  humid  air  in  the 
presence  CO  for  the  case  of  50  kV  are  shown  in  Figure  29.  To  analyze  the  concentration  profile  in 
this  case.  Equation  20  was  considered.  The  model  profile  obtained  in  this  case  is  very  similar  to  the 
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Figure  30  Comparison  of  the  Model  Prediction  of  NO  and  NO2  Profile  with  Ex¬ 
perimental  Data  in  Case  of  NO  Removal  in  the  Presence  of  500  ppm  of  CO  in  Humid 
Air 
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model  profiles  obtained  in  the  case  of  NO  removal  in  humid  air  (see  Figure  23).  The  addition  of  a 
single  reaction  of  carbon  monoxide  does  not  significantly  affect  the  species  concentration  profiles 
at  this  initial  CO  concentration. 

Figure  30  compares  the  NO  and  NO2  profiles  in  the  presence  of  CO  predicted  by  the  model 
to  the  experimental  points.  As  seen  in  the  figure,  the  profile  of  NO  matches  well  with  the 
experimental  points.  The  NO2  profile  of  the  model  qualitatively  matches  the  experimental  points 
but  the  values  are  significantly  different.  The  maximum  concentration  predicted  by  the  model  is 
around  50  ppm  whereas  the  maximum  observed  concentration  of  NO2  was  around  22  ppm.  The 
addition  of  high  concentrations  of  CO  (>  0.5%)  to  consider  the  removal  of  NO  may  be  considered 
for  future  work. 

G.  NO  REMOVAL  FROM  AIR  IN  THE  PRESENCE  OF  ETHYLENE 

The  atmospheric  chemistry  of  NO,^  and  hydrocarbons  has  been  well  characterized  by 
Atkinson  et  al.  (1984),  Seinfeld  (1986),  and  Grosjean  et  al.  (1996).  The  reactions  of  ethylene  in  the 
atmosphere  with  hydroxyl  radicals,  ozone,  and  dissociated  oxygen  may  also  be  applicable  in  the 
pulsed  streamer  corona  reactor. 

Mizuno  et  al.  (1993)  considered  the  effect  of  hydrocarbons  on  the  removal  efficiency  of  NO^ 
from  the  exhaust  of  a  diesel  engine.  They  observed  that  at  a  gas  temperature  of  30“C  the  NO^ 
removal  efficiency  was  22%  for  an  input  power  of  37.5  W.  With  the  addition  of  500  ppm  of 
ethylene,  a  rapid  increase  in  the  NO,(  removal  rate  (60%  for  an  input  power  of  35  W)  was  observed. 
Among  the  additives  considered,  ethylene  had  the  greatest  effect  on  the  removal  efficiency  of  NO,; 
compared  with  other  additives  such  as  methane  and  propane. 

Vogtlin  and  Penetrante  (1993)  conducted  experiments  on  the  removal  of  NO,;  in  the  presence 
of  hydrocarbon  additives  such  as  ethylene  and  octane  at  high  gas  temperatures  of  100  to  200“C. 
Octane  was  the  hydrocarbon  used  as  an  additive  at  a  concentration  ranging  from  0.8  to  8  times  the 
initial  NO,;  concentration.  At  a  ratio  of  8:1, 93.75%  of  the  initial  NO  was  removed  and  45.5%  of  the 
initial  NO,;  was  removed  at  an  energy  cost  of  29  eV  per  NO,;  molecule  removed.  When  the  molar 
ratio  was  lowered  to  4:1  for  a  gas  composition  of  5%  O2, 10%  H20, 15%  CO2,  and  70%  N2  (flue  gas 
composition),  almost  100%  NO  was  removed  and  45%  of  NO^  was  removed  at  an  energy 
consumption  of  80  eV  per  removed  NO,;  molecule.  No  byproducts  identification  was  made  during 


the  study.  The  hydrocarbon  additives  were  found  to  be  recycling  the  hydroxyl  radicals  dunng  the 
oxidation  and  reduction  of  NO.  Hence,  the  efficiency  of  a  particular  hydrocarbon  was  therefore 
limited  by  the  reaction  rate  of  this  hydrocarbon  with  hydroxyl  radicals. 

Wolf  et  al.  (1996)  also  considered  the  effects  of  hydrocarbons  on  NO,,  removal  m  a  dielectnc 
barrier  discharge  reactor.  They  found  that  the  addition  of  500  ppm,  1000  ppm,  and  2000  ppm 
ethylene  greatly  increased  the  removal  of  500  ppm  NO  from  humid  gas  containing  10%  water,  77 
%  nitrogen,  and  13  %  oxygen.  They  observed  only  a  small  amount  of  ethylene  decomposition  under 
these  conditions.  They  found  similar  effects  for  the  addition  of  toluene. 

In  the  natural  atmosphere,  ethylene  reacts  primarily  with  hydroxyl  radicals  since  the 
concentration  of  ozone  and  dissociated  oxygen  are  very  small.  Hydroxyl  radicals  directly  attack 
ethylene  to  produce  peroxyl  radicals  which  further  react  with  oxygen  molecules  through  the 
following  reactions: 

-  HOCHfH^-  (23) 

-  HOCHfHp^-  (24) 


NO  is  converted  to  NO,  through  the  following  reaction: 


HOCH^CHp^-*NO  -  NO^+HOCH^CHjO- 


(25) 


The  HOCHXH,0  radical  is  known  to  react  with  oxygen  to  give  glycolaldehyde  (HOCH^CHO) 
and  formaldehyde  (CH,0).  The  formaldehyde  produced  in  the  above  reactions  can  further  react 
with  hydroxyl  radicals  to  produce  hydroperoxy  radicals.  Direct  reactions  of  ethylene  with  ozone 
also  leads  to  the  production  of  formaldehyde,  formic  acid,  carbon  dioxide,  and  carbon  monoxide 
in  the  atmosphere  (Seinfeld,  1986)  tlirough  the  reaction 
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-  HCHO+OAHCOOH+O.ISCO^+OAICO 
O.UH^  +0A2HP+0. 1 IHO^' 


(26) 


Grosjean  et  al.  (1996)  studied  the  reactions  of  ozone  with  ethylene  and  found  that  under 
atmospheric  conditions  ethylene  reacts  with  ozone  to  produce  primarily  (>99  %)  formaldehyde 
with  very  trace  amounts  of  glyoxal  and  cyclohexanone. 

Ethylene  also  reacts  with  dissociated  oxygen  and  hydrogen  through  the  following 

reactions: 

C2/f4+0-  -  CH^-+HCO-  (27) 


C^H^+0-  -  CH^+CHp 


(28) 


C^H,+0-  -  C^H^-+OH- 


(29) 


(30) 


The  radicals  such  as  HCO-,  CHj-,  and  HOCH.CH,.  react  with  NO  and  NO,  to  form  other  products. 

DRY  AIR 

A  series  of  NO^  removal  experiments  at  different  voltages  (30,  35,  40,  and  50  kV)  and 
residence  times  (7.2, 14.5. 29.0,  and  44.0  seconds)  with  500  ppm  of  ethylene  present  in  dry  air  was 
considered. 

In  dry  air,  the  main  reaction  mechanism  for  ethylene  breakdown  was  through  reactions  with 
oxygen  radicals  to  produce  other  radicals  such  as  HCO-,  CH-,  CH2.  and  CH,0-.  Figure  3 1  shows  the 
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variation  of  NO  concentration  as  a  function  of  residence  time  for  various  operating  voltages.  As  can 
be  seen  in  the  figure,  the  NO  concentration  quickly  decreased  at  higher  voltages.  When  Figures  15 
and  3 1  are  compared,  it  can  be  seen  that  the  presence  of  ethylene  greatly  enhances  NO  removal.  NO 
removal  in  the  case  of  dry  air  and  in  the  absence  of  ethylene  at  7 .2  seconds  residence  time  and  at  an 
operating  voltage  of  50  kV  was  80%.  The  addition  of  ethylene  increases  the  removal  efficiency  to 
almost  100%  at  the  same  operating  conditions.  The  CH3.  and  HCO-  radicals  produced  from  ethylene 
breakdown  react  with  NO,  to  form  a  variety  of  compounds  such  as  CHjNO-  and  HNO-  thereby 
enhancing  the  NO  removal.  Also,  HCO-  radicals  increase  the  production  of  hydroperoxy  radicals, 
which  are  known  to  enhance  the  NO  removal: 


HCO-+0- 

-  //•+CO2 

(31) 

H‘+0^ 

-  HO^- 

(32) 

Figure  32  shows  the  NO,  concentration  profile  as  a  function  of  residence  time  for  various 
operating  voltages.  NO,  concentration  increased  monotonically  with  residence  time  for  the  lower 
voltages  (30  &  35  kV),  and  at  higher  voltages  (40  &  50  kV)  the  concentration  passes  through  a 
maximum.  At  the  highest  voltage  of  50  kV  and  at  longest  residence  time  of  44.0  seconds,  45% 
removal  of  NO,  is  observed  in  the  presence  of  500  ppm  of  ethylene  as  compared  to  5%  removal  of 
NO.  in  dry  air  without  ethylene.  When  Figures  16  and  32  are  compared  the  NO.  concentration  in 
the  presence  of  ethylene  (Figure  32)  is  significantly  smaller  at  higher  voltages  (40  &  50  kV)  and  at 
longer  residence  times  as  compared  to  the  profile  in  Figure  16  (no  ethylene).  The  HCO-  radicals 
react  with  NO,  to  form  HNQ  .  The  HCO-  radicals  also  produce  -OH  radicals  by  reaction  with 
oxygen  radicals.  The  corresponding  reactions  are: 

HCO-+NO^  -  HNO^^CO  (33) 
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Figure  31  NO  Concentration  Profile  in  an  Atmosphere  of  Dry  Air  and  500  ppm 
of  Ethylene 


Figure  32  NO2  Concentration  Profile  in  an  Atmosphere  of  Dry  Air  and  500  ppm 
of  Ethvlene 


The  hydroxyl  radicals  produced  by  Reaction  34  help  in  the  breakdown  of  ethylene  and  also 

HC0-+0-  -  CO+OH’  (34) 

in  the  removal  of  NO2.  Thus,  HCO-  is  a  very  important  reactive  radical  which  acts  as  a  pathway  for 
hydroxyl  and  hydroperoxy  radicals  and  also  acts  as  a  source  for  the  removal  of  NO,.  Some  of  the 
NO,  is  also  reduced  by  the  methyl  radicals  to  NO.  But.  the  main  reaction  mechanism  for  the  removal 

of  NO2  is  through  conversion  to  HNO,. 

A  list  of  additional  reactions  that  are  considered  for  modeling  NO  removal  in  the  presence 
of  ethylene  is  shown  in  Table  5.  The  model  concentration  profiles  for  various  chemical  species  for 
the  case  of  50  kV  are  shown  in  Figures  33  and  34.  Figure  33  shows  the  concentration  profiles  of 
NO,  NO,,  HNO,,  and  O3  as  functions  of  residence  time  along  with  the  experimental  data  of  NO  and 
NO,  for  the  case  of  50  kV  .  The  rate  constants  for  the  oxygen  and  nitrogen  dissociation  are  assumed 
to  be  the  same  as  those  calculated  from  the  NO  experimental  data  for  the  dry  air  case.  As  shown  in 
Figure  33  there  is  a  significant  formation  of  HNO2,  and  the  experimentally  observed  NO  removal 
was  faster  than  that  predicted  by  the  model.  At  a  residence  time  of  7.2  seconds  for  an  operating 
voltage  of  50  kV,  70%  removal  of  NO  is  predicted  by  the  model  as  compared  to  100%  removal  of 
NO  found  in  the  experiment.  In  addition,  the  model  shows  10  ppm  of  NO2  for  a  voltage  of  50  kV 
at  a  residence  time  of  29.0  seconds.  The  amount  of  NO,  coixesponding  to  the  experiment  at  the  same 
operating  conditions  is  91.6  ppm.  There  is  a  significant  difference  between  the  model  and  the 

experiments  in  this  case. 

Figure  34  shows  the  main  products  expected  firom  ethylene  breakdown  in  the  pulsed  streamer 
corona  reactor  that  were  predicted  by  the  model.  The  concentrations  of  CO  and  CH,0  increase 
Steadily  with  residence  time. 

HUMID  AIR 

A  series  of  experiments  at  different  voltages  (30,  35, 40,  and  50  kV)  and  residence  times  (7.2, 
14  5  29.0,  and  44.0  seconds)  was  conducted  to  consider  the  effect  of  ethylene  on  NO  removal  in 
humid  air.  The  presence  of  water  leads  to  the  corona-induced  production  of  hydroxyl  radicals  which 
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Mole  Fraction 


Figure  33  Concentration  Profile  of  Different  Species  Predicted  by  the  Model  for 
NO  Removal  in  the  Presence  of  500  ppm  Ethylene  in  Dry  Air  W  ith  the  Experimental 

Data 


Figure  34  .  Concentration  Profile  of  Breakdown  Products  Predicted  by  the  Model 
for  NO  Removal  in  Dry  Air  with  500  ppm  Ethylene 
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should  lead  to  more  ethylene  removal  than  that  observed  in  dry  air.  Figure  35  shows  the  profiles 
for  ethylene  breakdown  in  the  cases  of  dry  and  humid  air  as  predicted  by  the  model. 

The  concentration  profiles  of  NO  as  a  flmction  of  residence  time  in  humid  air  with  ethylene 
for  various  operating  voltages  are  shown  in  Figure  36.  The  rate  of  NO  removal  is  lower  compared 
to  the  profile  obtained  for  dry  air  with  ethylene  (Figure  3 1).  The  concentration  of  NO  at  7.2  seconds 
residence  time  was  17  ppm  in  humid  air  as  opposed  to  almost  0  ppm  m  the  dry  air  (both  with 

ethylene). 

The  concentration  profiles  of  NO^  in  humid  air  with  ethylene  are  shown  in  Figure  37.  The 
NO,  concentration  profiles  go  through  maxima  similar  to  the  NO,  profiles  in  humid  air  without  the 
presence  of  ethylene  (Figure  25).  The  maximum  concentration  of  NO,  at  50  kV  was  about  24.4  ppm 
at  a  gas  residence  time  of  7.2  seconds,  which  compares  well  with  the  NO,  maximum  concentration 
of  21.2  ppm  in  humid  air  without  ethylene.  Hence,  ethylene  does  not  significantly  affect  NO, 
removal  in  the  presence  of  water.  The  model  profiles  indicate  that  the  main  breakdown  product 
formed  from  NO,  in  the  case  of  humid  air  and  500  ppm  ethylene  is  HNO,  as  opposed  to  HNO3  in 
humid  air  with  no  ethylene.  More  extensive  experimental  analysis  of  nitrite  and  nitrate  products  is 

needed  to  confirm  this  prediction. 

The  concentration  profiles  of  NO,  NO,,  HNO,.  and  HNQ  in  humid  air  with  500  ppm  of 
ethylene  predicted  by  the  model  are  shown  in  Figure  38  along  with  the  experimental  data  of  NO  and 
NO,.  The  model  profile  of  NO  matches  very  well  with  the  experimental  data.  The  profile  of  NO, 
predicted  by  the  model  also  qualitatively  matches  with  the  experimental  data.  The  model 
concentration  profiles  of  the  main  byproducts  formed  from  ethylene  breakdown  such  as 
glycolaldehyde  (HOCH,CHO),  formaldehyde  (CH,0),  and  CO  are  shown  in  Figure  39.  To  calculate 
the  profile,  the  rate  constants  for  nitrogen,  oxygen,  and  water  dissociation  have  been  assumed  to  be 

the  same  as  that  obtained  for  NO  removal  from  humid  air. 

To  confirm  the  validity  of  the  model,  a  series  of  experiments  to  determine  the  byproducts 
formed  by  the  dissociation  of  ethylene  in  dry  air  was  conducted  at  50  kV  and  at  44.0  seconds 
residence  time.  The  gas  sample  analysis  by  GC/MS  indicated  the  complete  breakdown  of  ethylene. 
For  500  ppm  ethylene  one  would  expect  1000  ppm  CO,  if  complete  conversion  were  achieved; 
however,  130  ppm  of  CO,  was  observed  thereby  indicating  only  partial  degradation.  HPLC  analysis 
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Figure  35  ;  Ethylene  Model  Profile  With  and  Without  the  Presence  of  Water  in 
C^se  of  NO  Removal  from  Air 


Figure  36  :  NO  Concentration  Profile  in  an  Atmosphere  of  Humid  Air  and  500  ppm 
of  Ethylene 
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Figure  37  ;  Concentration  Profile  of  Breakdown  Products  Predicted  by  the  Model 
in  Case  of  NO  Removal  in  Humid  Air  with  500  ppm  of  Ethvlene 


Table  4  '  Byproducts  Predicted  by  the  Model  Due  to  Ethylene  Breakdown 


Name 

Chemical  Formula 

Formaldehyde 

CHoO 

Formic  Acid 

HCOOH 

Glycolaldehyde 

HOCH.2CHO 

Acetaldehyde 

CH3CHO 

Ethylene  Glycol 

OHCH2CH2OH 

Carbon  Monoxide 

CO 

Carbon  dioxide 

<M 

0 

0 

Methyl  Nitrate 

CH3NO 

Hydrogen  Cyanide 

HCN 
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of  the  liquid  samples  from  dry  air  showed  no  peaks,  however,  analysis  of  humid  air  reactor  samples 
showed  four  peaks  (Figure  40).  One  of  these  peaks  was  positively  identified  to  be  formic  acid. 
Mizuno  et  al.  (1995)  considered  the  removal  of  NO,;  gases  from  dry  air  and  found  acetic  acid  to  be 
an  important  byproduct.  Acetic  acid  was  not  identified  as  one  of  the  byproducts  m  the  present 
analysis,  nor  is  it  predicted  by  any  of  the  reaction  mechanisms  from  the  atmospheric  chemistry 
literature.  Table  5  shows  the  list  of  chemical  compounds  that  could  be  formed  from  the  ethylene 

removal. 
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Figure  38  NO,  Concentration  Profiles  in  an  Atmosphere  of  Humid  Air 
and  500  ppm  Ethylene. 


Figure  39  Concentration  Profiles  of  Various  Species  Predicted  by  the 

Model  in  the  Case  of  NO  Removal  in  the  Presence  of  500  ppm  Ethylene 
in  Humid  Air. 
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Figure  40.  ChroTnatogram  from  HPLC  Analysis  of  Gas 
from  Corona  Experiment  at  50  kV  with  NO,  ethylene , 
air 
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Table  5  :  Additional  List  of  Chemical  Reactions  Considered  for  Modeling  NO  Re¬ 
moval  in  the  Presence  of  Ethylene 


Chemical  Reactions 

Rate  Constants  Source 

{cm^  .mole~^  .sec~^) 

C2H4  +  0--^CH,  +  HC0- 

4.39  X  10^^ 

NIST 

C2H4  +  0-  — ^  CH2  +  CHiO 

3.98  X  10^^ 

NIST 

Gift  +0-^H+  CH,CHO 

3.76  X  10^^ 

NIST 

C2H4  +  0-  — >  C2H3  +  OH- 

1.46  X  10^ 

NIST 

C2H4  +  N-  — ^  HCX  +  CH2 

9.93  X  10^° 

NIST 

C2H4  +  H-  — >  C2H-^ 

6.88  X  10^^ 

NIST 

C2H4  +  OH-  — ^  C2H4OH- 

5.42  X  10^2 

NIST 

HCO-  +  CH^  — ^  CO  ^  CH4 

1.21  X  10^^ 

NIST 

HC0-  +  0-  — ^  CO  +  OH- 

3.10  X  10^^ 

NIST 

HCO-  +  0-  — ^  H-  +  CO2 

3.01  X  10^2 

NIST 

HCO-  OH-  — ^  H2O  +  CO 

1.02  X  10^“^ 

NIST 

HCO-  +  O2  — ^  HO-.,  +  0- 

3.31  X  10^2 

NIST 

HCO-  +  NO2  — ^  CO  +  HNO2 

3.27  X  10^^ 

NIST 

HCO-  +  NO  — ^  CO  d-  HNO 

7.23  X  10^2 

NIST 

HCO-  -t-  Ci?3  CH3CHO 

1.81  X  10^^ 

NIST 

C2H4OH-  +  O2  — ^  HOCH2CH2O-2 

1.81  X  10^2 

NIST 

2HOCH2CH2O-2 

2HOCH2CH2O-  +  O2 

5.0  X  10^^ 

NIST 

continued  on  next  page 
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continued  from  previous  page _ _ 

Chemical  Reactions 

2HOCH2CH2O., 

HOCH2CH2OH  +  HOCHoCHO  +  O2 
H0CH2CH20-^  +  NO 
HOCH2CH20-  +  NO2 
CHi  +  HO2  — >  CHzO-  +  OH 
CH,  +  0—  CH^O 
CH^  +  0-  H-  +  CH2O 
CH2  +  H0:2  — >  CH4  +  O2 
CH,  +  A'O  —  CH2N0 
CH^  +  NO2  — ^  CH3ONO 
CH^  +  NO2  — ^  CH^O-  +  NO 
CH20  +  0-  — ^  OH- +  HCO- 
CHiO  +  OH-  — ^  H2O  +  HCO- 
CH2O  +  OH-  — >  HCOOH  ^  H- 
CH2  +  OH-  CH2O  +  H- 
CH2  +  HCO-  — ^  CHi  +  CO 
CHiCHO  +  O2  — ^  CH2O  ^CO-irOH- 
C2H,  +  0-  -^CHsCHO  +  H- 
C2H-3  +  HCO-  — ^  C2H-,CH0 

C2i?5  +  HCO-  - >•  ^2-^6  +  CO 

+  OH-  CH3CHO 
CH^O-  +  HCO-  — ^  CH3OH  +  CO 


Rate  Constants  Source 
[cm^  .mole~^  .sec~^) _ 


4.04  X  10^^ 

NIST 

5.42  X  10^^ 

NIST 

1.81  X  10^^ 

NIST 

1.58  X  10^° 

NIST 

8.4.3  X  10^' 

NIST 

3.61  X  10^^ 

NIST 

1.01  X  10^2 

NIST 

7.00  X  10^° 

NIST 

1.39  X  10^'^ 

NIST 

1.01  X  10" 

NIST 

6.03  X  10" 

NIST 

1.21  X  10" 

NIST 

1.81  X  10" 

NIST 

1.81  X  10" 

NIST 

1.81  X  10" 

NIST 

8.02  X  10" 

NIST 

1.81  X  10" 

NIST 

1.21  X  10" 

NIST 

3.01  X  10" 

NIST 

9.04  X  10" 

NIST 

continued  on  next  page 
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continued  from  previous  page 

Chemical  Reactions 

Rate  Constants 

(cm^  .mole~^  .sec~^) 

Source 

E-  +  CH3OH  — ^  CH2OH  +  H2 

8.39  X  10^2 

NIST 

0-  +  CHzOH  — ^  OH-  +  CH2OH 

4.87  X  10® 

NIST 

OH-  +  CH3OH  CH2OH  +  H2O 

6.4  X  10“ 

NIST 

CHzO-  +  CH3OH  — ^  CHzOH  +  CH2OH 

3.1  X  10® 

NIST 

NO3  +  CH3OH  — >  HNO3  +  CH2OH 

5.73  X  10® 

NIST 

CH2OH  +  CH2OH  HOCH2CH2OH 

4.82  X  10“ 

NIST 
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SECTION  V 
CONCLUSIONS 


In  the  present  study,  it  was  shown  that  a  pulsed  streamer  corona  discharge  can  be  used  to 
remove  NO^  from  feed  gases  such  as  nitrogen  and  air  containing  small  amounts  of  other  gases  such 
as  water,  carbon  monoxide,  and  ethylene.  A  rotating  spark  gap  pulsed  power  supply  providing  fast 
rise  time,  narrow  width,  repetitive  high  voltage  pulses  to  a  wire-cylinder  geometry  gas  phase  reactor 
vessel  was  constructed  and  employed  in  this  study.  Experiments  varying  the  gas  composition,  the 
pulsed  voltage  value,  and  the  gas  residence  time  were  performed  to  determine  the  NO  and  NO2 
removal.  The  main  reaction  pathways  for  NO  removal  were  identified  for  each  of  the  gas 
compositions.  A  model  was  developed  to  calculate  the  concentration  profiles  of  different  chemical 
species  under  treatment  by  a  pulsed  streamer  corona  discharge.  The  reaction  rate  constants  for  the 
nitroaen.  oxygen,  and  water  molecule  dissociations  were  obtained  by  fitting  the  experimental  data 
of  NO  removal  with  the  model. 

Power  input  into  the  pulsed  corona  discharge  reactor  was  calculated  using  the  voltage  and 
current  measurements  made  during  experiments  with  a  fast  oscilloscope  and  high  voltage  and  current 
probes  attached  to  the  pulsed  power  input.  The  measured  peak  voltage  was  significantly  higher  than 
the  applied  voltage,  and  the  pulse  rise  time  varied  from  15-30  nanoseconds  as  the  applied  voltage 
was  increased  from  30  to  50  kV.  The  pulse  width  at  half  maximum  varied  from  500  to  700 
nanoseconds.  The  energy  input  varied  from  8.0  mJ/pulse  at  30  kV  to  148  mJ/pulse  at  50  kV  in  an 
atmosphere  of  dry  air.  The  power  input  varied  from  0.48  W  to  8.88  W  at  the  corresponding 

conditions. 

Figure  41  shows  the  summary  of  NO  concentration  profiles  in  different  atmospheres  at  a 
voltage  of  50  kV,  and  Figure  42  shows  the  summary  of  the  corresponding  NOj  concentration 
profiles.  As  seen  in  Figure  41,  NO  removal  in  the  presence  of  ethylene  in  dry  air  was  the  fastest. 
The  addition  of  water  did  not  enhance  the  NO  removal.  The  concentration  profiles  of  NO  removal 
in  humid  air,  in  the  presence  of  500  ppm  of  CO  in  humid  air,  and  in  the  presence  of  500  ppm  of 
ethylene  in  humid  air  almost  match  each  other.  The  NO.  concentration  profiles  shown  in  Figure  42 
indicate  that  in  the  presence  of  water,  NO.  formation  is  significantly  reduced  compared  to  NQ 
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removal  in  dry  air.  The  addition  of  ethylene  to  dry  air  enhanced  the  NO,  removal  as  shown  in 
Figure  42.  As  was  observed  in  the  case  of  NO  removal,  the  NO2  concentration  profiles  in  the  case 
of  humid  air,  500  ppm  of  CO  in  humid  air,  500  ppm  of  ethylene  in  humid  air  matched  each  other 

well. 

For  the  experimental  case  of  pure  nitrogen  containing  NO  under  pulsed  corona  treatment, 
NO  was  chemically  reduced  to  N^.  About  99%  removal  of  NO,  was  observed  at  50  kV  pulsed 
voltage  for  a  gas  residence  time  of  44.0  seconds.  The  model  profiles  of  NO  and  NOj  matched  well 
with  the  experimental  data,  and  the  model  indicated  that  there  was  a  small  amount  of  N2O  formation 

during  the  process. 

In  the  experimental  case  of  dry  air  containing  NO,  oxidation  mechanisms  dominate  the 
reduction  pathway  and  hence  NO  was  converted  to  NO,.  Subsequent  removal  of  NO,  was  not 
observed  under  the  current  operating  conditions  (maximum  of  50  kV  pulsed  voltage  and  60.0 
seconds  residence  time).  The  addition  of  oxygen  in  dry  air  as  the  feed  gas  increased  the  rate  of  NO 
removal.  Chang  et  al.  (1992)  reported  that  in  a  dielectric  barrier  discharge  system,  the  addition  of 
oxygen  decreased  the  rate  of  NO  removal.  The  model  developed  in  the  present  study  predicted  80% 
conversion  of  the  initial  NO  to  NO^,  while  the  experimental  value  showed  95%  conversion  of  NO 
to  NO2  at  50  kV  and  at  long  residence  time. 

The  rate  constant  for  nitrogen  dissociation  was  determined  by  fitting  the  model  with  the  NO 
experimental  data.  The  experiments,  as  well  as  the  modeling  results,  match  well  with  the  results 
indicated  in  the  non-thermal  plasma  treatment  literature.  The  rate  constants  for  nitrogen  and  oxygen 
dissociation  were  determined  from  the  NO  removal  and  the  NO2  formation  data.  The  reaction  rate 
constant  for  nitrogen  dissociation  was  found  to  be  different  from  the  value  obtained  in  pure  nitrogen. 
The  rate  constant  was  0.0475  s  ''  in  an  atmosphere  of  nitrogen  at  50  kV,  but  in  dry  air  the  rate 
constant  was  0.03  s  Hence,  as  expected,  the  rate  constant  depends  upon  the  composition  of  the 

gas. 

The  N2O  measurement  matched  well  with  the  model  value  obtained  in  the  experiment.  This 
result  is  significant  since  this  is  an  independent  check  of  the  model.  To  have  an  additional 
independent  check  on  the  model,  experiments  were  conducted  to  measure  ozone  formation  in  dry 
air  under  pulsed  streamer  corona  discharge.  The  model  prediction  matched  qualitatively  with  the 


79 


experimental  value  for  ozone,  however,  further  work  is  necessary  to  fully  validate  the  experimental 
ozone  measurements. 

In  the  case  of  NO  removal  from  feed  air  in  the  presence  of  water,  the  addition  of  water  did 
not  enhance  NO  removal.  The  NO2  concentration  was  significantly  lower  due  to  the  water-induced 
formation  of  further  products  from  NO2  reactions.  A  very  small  amount  of  HNO3  was  measured  as 
nitrate  during  the  experiments  using  ion  chromatography.  The  model  predicts  about  35%  conversion 
of  the  initial  NO  to  HNO3  at  50  kV  and  at  long  residence  time  (60.0  seconds).  The  experimental 
measurements  did  not  find  such  high  concentrations  of  nitric  acid.  The  NO2  concentration  profile 
qualitatively  matched  with  the  experimental  data,  but  quantitatively  the  difference  was  significant. 

The  presence  of  500  ppm  of  carbon  monoxide  in  humid  air  did  not  affect  the  experimental 
concentration  profiles  of  NO  and  NOj.  Tokunaga  et  al.  (1984)  showed  that  the  removal  of  NO  was 
enhanced  in  the  presence  of  1.13%  carbon  monoxide  by  volume.  In  the  present  study  the  amoimt 
of  carbon  monoxide  added  to  the  gas  may  have  been  too  small  to  have  an  affect  on  the  concentration 
of  NO.  The  model  profile  with  and  without  500  ppm  of  carbon  monoxide  also  showed  no  effect  of 
CO  on  NO  removal. 

In  the  presence  of  500  ppm  of  ethylene,  NO  removal  was  significantly  enhanced.  In  dry  air 
with  ethylene  about  45%  removal  of  NO2  was  observed.  The  main  reaction  pathways  for  the 
enhanced  removal  of  NO  and  NO2  are  due  to  the  production  of  highly  reactive  radicals  by  the 
reaction  of  oxygen  atoms  with  ethylene.  The  formation  of  radicals  such  as  HCO-,  and  CH3-  react 
with  NO  and  NO2  to  enhance  their  removal.  The  model  indicated  that  the  main  reaction  pathway 
for  the  removal  of  NO2  in  dry  air  is  through  reaction  with  HCO-  radical  to  form  HNO2. 

In  humid  air  containing  ethylene,  additional  ethylene  breakdovm  is  caused  by  its  reaction 
with  hydroxyl  radicals,  aiding  in  the  removal  of  NO2.  Complete  removal  of  ethylene  under  these 
conditions  was  found  by  using  GC/MS  analysis  of  the  gas  samples.  Four  chromatographic  peaks 
were  observed  when  the  sample  was  analyzed  by  HPLC.  Formic  acid  was  identified  as  one  of  these 
peaks.  The  amount  of  formic  acid  represents  less  than  1  %  of  the  carbon  initially  fed  into  the  reactor. 
The  production  of  carbon  dioxide  was  observed  for  dry  air  conditions;  however,  quantitatively  only 
13%  of  the  input  carbon  could  be  accounted  for  by  CO2. 

In  conclusion,  the  use  of  a  pulsed  streamer  corona  discharge  has  been  shown  to  be  successful 
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in  the  removal  of  NO^  in  different  gas  compositions.  The  main  reaction  mechanisms  for  1^0 
removal  in  various  gas  compositions  under  corona  discharge  were  identified  in  this  work.  The 
addition  of  ethylene  was  found  to  enhance  both  NO  and  NO2  removal  even  in  the  case  of  dry  air, 
while  water  addition  assisted  NOj  but  not  NO  removal.  The  addition  of  500  ppm  of  CO  in  humid 
air  did  not  enhance  NO  and  NO2  removal. 
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Figure  41  :  Summary  of  NO  Removal  in  Different  Atmospheres 


Figure  42;  Summary  of  NO2  removal  in  Different  Atmospheres 
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SECTION  VI 
RECOMMENDATIONS 


Recommendations  from  this  study  include: 

1)  Aerosol  Measurements.  Experiments  conducted  to  study  the  removal  of  NO^  in  humid  air 
indicated  both  NO  and  NOj  removal.  The  atmospheric  chemistry  literature  indicates  that  both 
HNOj  and  HNOj  would  be  produced  and  could  form  aerosols.  The  analysis  of  aerosols  is  an 
important  area  for  future  work  in  order  to  develop  a  more  complete  material  balance  on  the 
products  of  reaction.  Methods  to  characterize  the  size  distribution  of  the  aerosols  and  to  measure 
the  concentration  of  aerosols  would  be  needed  for  these  studies. 

2)  Operating  Conditions.  Variation  of  the  gas  temperature,  and  additional  gas  compositions 
should  be  considered.  The  exhaust  gases  emitted  from  power  plants  and  jet  cell  engines  are  at 
high  temperature  and  contain  large  concentrations  of  carbon  dioxide,  water  vapor,  and  aerosol 
particles.  The  addition  of  aerosols  and  operation  of  the  reactor  at  higher  temperatures  may  be  of 
interest.  Further  variation  of  the  gas  composition  through  consideration  of  higher  concentrations 
of  NO,  CO,  CO2,  water  vapor,  and  other  could  be  considered. 

3)  Improvements  in  Experimental  Setup.  The  power  input  into  the  system  is  currently  adjusted 
by  increasing  the  applied  voltage.  Under  the  current  operating  conditions,  the  maximum  voltage 
achievable  is  50  kV.  An  alternative  method  to  increase  the  power  input  into  the  system  would  be 
to  vary  the  frequency  of  pulsing.  Modification  of  the  power  supply  to  allow  for  frequency 
variation  would  be  useful  in  order  to  further  evaluate  the  difference  between  the  effects  of 
applied  field  strength  and  power  density  on  NO  removal. 

4)  Reactor  electrode  materials.  It  has  been  observed  in  the  literature  that  in  order  to  achieve 
reducing  conditions  in  the  reactor,  it  is  desirable  to  run  the  reactor  with  an  internal  source  of 
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carbon.  It  may  therefore  be  of  interest  to  operate  the  reactor  with  a  packed  bed  of  activated 
carbon  or  other  carbon  particles.  In  addition,  the  use  of  carbon  electrodes  may  also  be  of  interest 
in  order  to  facilitate  degradation  reactions. 

5)  Improvements  in  Modeling.  It  is  currently  assumed  in  the  modeling  work  that  the  electric 
field  is  constant  during  the  pulse-on  period.  Variation  of  the  electric  field  pulse  shape  during  the 
pulse-on  period  and  variation  along  the  electrode  length  could  be  introduced  to  acount  for  non- 
uniform  aspects  of  the  electric  discharge. 
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APPENDIX  A 


NO  REMOVAL  DATA  FROM  NITROGEN,  DRY,  AND  HUMID  AIR 


The  concentrations  shown  in  the  tables  have  been  normalized  for  an  initial  con¬ 
centration  of  NO  to  be  100  ppm  and  NO2  concentration  to  be  0  ppm. 

Table  A.l:  NO  and  NO2  Data  at  30  kV  in  an  Atmosphere  of  Nitrogen 


Res.  Time 

(sec) 

NO  Cone. 

(ppm) 

NO2  Cone. 

(ppm) 

Mean  NO 

Mean  NO2 

Runl 

Run2 

Runl 

Run  2 

7.2 

99.8 

98.8 

0 

0 

99.9 

0 

14.5 

87.6 

87.7 

3.3 

4.8 

87.6 

4.0 

29.0 

82.7 

80.7 

4.1 

4.5 

81.7 

4.3 

44.0 

57.9 

64.5 

7.0 

7.9 

61.2 

7.4 
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Table  A. 2:  NO  and  NO2  Data  at  35  kV  in  an  Atmosphere  of  Nitrogen 


Res.  Time 

NO  Cone. 

NO2 

Cone. 

Mean  NO 

Mean  NO2 

(sec) 

(ppm) 

(ppm) 

Runl 

Run  2 

Runl 

Run  2 

7.26 

98.46 

90.74 

2.6 

3.0 

94.62 

2.80 

14.52 

64.8 

65.43 

6.94 

7.2 

65.11 

7.06 

29.04 

43.85 

44.96 

8.38 

7.02 

44.05 

4.93 

44 

28.75 

28.55 

6.76 

7.36 

28.65 

6.56 

Table  A. 3:  NO  and  NO2  Data  at  40  kV  in  an  Atmosphere  of  Nitrogen 


Res.  Time 

NO  Cone. 

NO2 

Cone. 

Mean  NO 

Mean  NO2 

(.sec) 

(ppm) 

(ppm) 

Runl 

Run2 

Runl 

Run  2 

7.2 

82.75 

78.05 

5.17 

7.03 

80.42 

6.10 

14.5 

■ 

47.29 

46.95 

7.78 

8.06 

47.12 

7.92 

8.96 

4.46 

5.41 

10.91 

4.93 

44 

1.58 

1.98 

2.0 

2.06 

1.78 

2.30 

Table  A.4:  NO  and  NOj  Data  at  50  kV  in  an  Atmosphere  ot  Nitrogen 


Res.  Time 

NO  Cone. 

NO2  Cone. 

Mean  NO 

Mean  NO2 

(sec) 

(ppm) 

(ppm) 

Runl 

Run2 

Runl 

Run  2 

7.2 

48.07 

49.17 

7.78 

8.0  1 

48.62  ' 

7.87 

14.5 

12.36 

12.4 

2.8 

3.0 

12.38 

2.9 

29 

0 

0 

0.76 

0.6 

0 

0.68 

44 

0 

0 

0.83 

1.24 

0 

1.03 

Table  A,5;  NO  and  SOi  Data  at  30  kV  in  an  Atmosphere  of  Dry  Air 


Res.  Time 

NO  Cone. 

NO2  Cone. 

Mean  NO 

Mean  NO2 

(sec) 

(ppm) 

(ppm) 

Runl 

Run2 

Runl 

Run  2 

3.6 

97.5 

98.1 

0 

97.8 

0 

7.2 

97.2 

97.6 

0.65 

1.62 

97.4 

1.13 

14.5 

97.4 

99.2 

3.3^ 

0 

98.3 

1.6 

29.0 

79.8 

94.6 

28.0 

3.7 

87.2 

15.8 

59.0 

77.0 

73.0 

27.3 

22.8 

75 

25.0 

04/05/2000  12:43 
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Table  A.6:  NO  and  NOj  Dlta  35  kV  in  an  Atmosphere  of  Dry  Air 


Res,  Time 

(sec) 

NO  Gone. 

(ppm) 

NO2  Gone. 

(ppm) 

Mean  NO 

Mean  NO^ 

Runl 

Run^;' 

Runl 

Run  2 

99.0 

98.6 

0 

0 

98.7 

0 

7.2 

82.04 

87.1# 

17.5 

11.5 

84.58 

M..5 

14.5 

90.5 

95.4  ^ 

12.4 

6,4 

92.9 

9,4 

29 

47.4 

28.6|| 

03.4 

84.8 

38 

74.1 

59 

9.5 

106.1 

107.2 

10.1 

106.1 

Table  A. 7:  NO  and  NO*  Dita  at  40  kV  in  an  Atmosphere  of  Dry  Air 


Res.  Time 

(sec) 

NO  Cone. 

(ppm)  ^ 

NO3 

(pi 

Cone. 

pm) 

Mean  NO 

Mean  NOj 

Runl 

15 

Runl 

Run  2 

3.6 

98.9 

in 

m 

3.0 

97.1 

1.5 

58.56 

68.9ip 

42.1 

33.0 

0.3.8 

.17.0 

j■E9i 

44.4 

mm 

64.1 

64.7 

41.8 

()4.4 

29 

0 

mi 

121.9 

117.9 

1.9 

1 19.9 

1.1 

mi 

99.4 

101.9 

0.5 

100.6 

Table  A.8:  NO  and  NO2  Data  at  50  kV  in  an  Atmosphere  of  Dry  Air 


Res.  Time 

(sec) 

NO  Cone. 

(ppm) 

NO2  Cone. 

(ppm) 

Mean  NO 

Mean  NO2 

Runl 

Run2 

Runl 

Run  2 

3.6 

75.0 

95.7 

24.7 

3.3 

85.3 

14.0 

7.2 

18.8 

21.3 

82.5 

79.8 

20.0 

81.2 

14.5 

0 

1.9 

108.2 

108.6 

0.9 

108.4 

29 

0 

0 

111.9 

119.0 

0 

115.4 

44 

0 

0 

85.4 

87.4 

0 

86.4 

Table  A. 9:  NO  and  NO2  Data  at  30  kV  in  an  Atmosphere  of  Humid  Air 


Res.  Time 

(sec) 

NO  Cone. 

(ppm) 

NO2  Cone. 

(ppm) 

Mean  NO 

Mean  NO2 

Runl 

Run2 

Runl 

Run  2 

7.2 

97.64 

97.94 

1.0 

0 

97.79 

0.5 

14.5 

98.66 

94.0 

2.04 

3.62 

96.33 

2.83 

29.0 

80.64 

88.67 

2.42 

2.4 

84.66 

2.4 

44 

89.59 

91.82 

1.02 

0 

90.7 

0.5 

04/05/2000  12:43 
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Table  A, 10:  NO  and  NO2  at35  kV  in  an  Atmosphere  of  lluinid  Air 


lable  A. 11:  NO  and  NOj  Dj||[a:at  40  kV  in  an  Atmosphere  of  Humid  Air 


Table  A. 12:  NO  and  NO2  Data  at  50  kV  in  an  Atmosphere  of  Humid  Air 


Res.  Time 

(sec) 

NO  Cone. 

(ppm) 

NO2  Cone. 

(ppm) 

Mean  NO 

Mean  NO 2 

Runl 

Run2 

Runl 

Run  2 

7.2 

28.1 

22.06 

20.8 

21.58 

25.08 

21.19 

14.5 

0 

0 

4.75 

10.26 

0 

7.5 

29.0 

0 

0 

2.76 

2.62 

0 

2.69 

44.0 

0 

0 

3.81 

3.83 

0 

3.82 

APPENDIX  B 


NO  REMOVAL  DATA  WITH  THE  ADDITION  OF  CO  AND  C2H4 


The  concentrations  shown  in  the  table  have  been  normalized  for  an  initial  con 
centration  of  NO  to  be  100  ppm  and  NO2  concentration  to  be  0  ppm. 

Table  B.l:  NO  and  NO2  Data  at  30  kV  in  an  Atmosphere  of  Humid  Air  and  CO 


Res.  Time 

NO  Cone. 

NO2 

Cone. 

Mean  NO 

Mean  NO2 

(sec) 

(ppm) 

(ppm) 

Runl 

Run2 

Runl 

Run  2 

7.2 

98.6 

97.5 

0.2 

0 

98.0 

0.1 

14.5 

96.3 

88.8 

1.5 

4.0 

92.6 

2.8 

29.0 

94.2 

93.3 

1.2 

1.1 

93.8 

1.2 

44.0 

93.2 

92.6 

1.6 

0 

92.9 

0.8 

Table  B.2:  NO  and  NO2  Data  at  35  kV  in  an  Atmosphere  of  Humid  Air  and  CO 


Res.  Time 

(sec) 

NO  Cone. 

(ppm) 

NO2  Cone. 

(ppm) 

Mean  NO 

Mean  NO2 

Runl 

Run2 

Runl 

Run  2 

7.2 

94.6 

94.3 

1.6 

1.2 

94.5 

1.4 

14.5 

78.0 

60.8 

3.4 

9.0 

69.4 

6.2 

29.0 

50.2 

41.3 

6.4 

6.7 

45.8 

6.6 

44.0 

45.6 

42.6 

5.7 

6.1 

44.1 

5.9 

Table  B.3:  NO  and  NO2  Data  at  40  kV  in  an  Atmosphere  of  Humid  Air  and  CO 


Res.  Time 

(sec) 

NO  Cone. 

(ppm) 

NO2  Cone. 

(ppm) 

Mean  NO 

Mean  NO2 

Runl 

Run2 

Runl 

Run  2 

7.2 

66.3 

58.9 

13.6 

19.5 

62.6 

16.6 

14.5 

27.5 

26.0 

11.5 

12.3 

26.8 

11.9 

29.0 

5.6 

4.2 

7.1 

6.75 

4.9 

6.9 

44.0 

9.6 

2.5 

5.6 

5.9 

6.05 

5.8 

Table  B.4:  NO  and  NO2  Data  at  50  kV  in  an 


Atmosphere  of  Humid  Air  and  CO 


Res.  Time 

(sec) 

NO  Cone. 

(ppm) 

NO2  Cone. 

(ppm) 

Mean  NO 

Mean  NO2 

Runl 

Run2 

Runl 

Run  2 

7.2 

30.5 

28.4 

25.6 

24.3 

29.5 

24.9 

14.5 

4.6 

2.3 

10.8 

9.7 

3.5 

10.3 

29.0 

0 

0 

2.5 

2.9 

0 

2.7 

44.0 

0 

0 

5.2 

5.0 

0 

5.1 

Table  B.5:  NO  and  NO2  Data  at  30  kV  in  an  Atmosphere  of  Dry  Air  and  500  ppm 
Ethylene 


Res.  Time 

(sec) 

NO  Cone. 

(ppm) 

NO2  Cone. 

(ppm) 

Mean  NO 

Mean  NO2 

Runl 

Run2 

Runl 

Run  2 

7.2 

91.8 

97.9 

1.6 

0 

94.9 

0.8 

14.5 

92.4 

88.9 

7.6 

10.8 

90.65 

9.2 

29.0 

79.0 

70.4 

20.0 

20.1 

74.7 

20.1 

44.0 

27.9 

48.0 

68.5 

43.4 

.38.0 

.56.0 

Table  B.6:  NO  and  NO2  Data  at  35  kV  in  an  Atmosphere  of  Dry  Air  and  500  ppm 
of  Ethylene 


Res.  Time 

NO  Cone. 

NO2 

Cone. 

Mean  NO 

Mean  NO2 

(sec) 

(ppm) 

(ppm) 

Runl 

Run2 

Runl 

Run  2 

7.2 

65.1 

71.9 

36.8 

27.9 

68.5 

31.9 

14.5 

66.4 

77.9 

34.3 

23.9 

72.2 

29.1 

29.0 

42.6 

24.9 

54.0 

63.1 

33.8 

58.5 

44.0 

1 

1.1 

91.7 

87.0 

1.1 

89.4 

Table  B./:  NO  and  NO2  Data  at  40  kV  in  an  Atmosphere  of  Dry  Air  and  500  ppm 
Ethylene 


Res.  Time 

NO  Cone. 

NO2 

Cone. 

Mean  NO 

Mean  NO2 

(sec) 

(ppm) 

(ppm) 

Runl 

Run2 

Runl 

Run  2 

7.2 

25.5 

14.7 

77.2 

87.0 

20.1 

82.1 

14.5 

4.3 

3.7 

95.6 

94.8 

4.0 

95.2 

29.0 

0 

0 

90.0 

93.2 

0 

91.6 

44.0 

0 

0 

63.9 

57.3 

0 

60.6 

04/05/2000  12:43 
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r<il>lf:»  IJ,8;  NO  and  NO2  Data  at  60  kV  in  an  Atmoflphf re  of  Dry  Air  uikI  ■>()() 
I’lUiylnn' 


R.rs.  Time 

NO  Cone. 

NO2  Cone. 

Mean  NO 

Mean  NO-.. 

(see) 

(ppm) 

(ppm) 

Runl 

Run2 

Runl 

Run  2 

0 

0 

102.3 

99.0 

0 

100.7 

l't.5 

0 

0 

86.1 

84.0 

0 

8.5.1 

2fl.O 

0 

0 

88.6 

88.4 

0 

88-5 

4'1.0 

0 

0 

55.0 

59.5 

0 

57..’{ 

[able  liJ):  NO  and  NO2  Data  at  30  kV  in  an  Aliiiosphcrc  of  Himiid  Air  und  500 
p|)in  KUiyU^ip 


Res.  'rime 

NO  Cone. 

NOj 

Cone. 

Meati  NO 

Mean  NO2 

(sec) 

(ppm) 

(ppm) 

Runl 

Ruii2 

Rtml 

Run  2 

97.6; 

— 

0 

98.0 

0 

14.5 

76.1 

73.2 

8.3 

■W 

74.6.5 

8.6 

29.0 

92.2 

73.9 

1.0 

83.1 

2.;) 

H.O 

80.0 

86.8 

0 

0 

82.!) 

0 

Table  B.IO:  NO  and  NO2  Data  at  35  kV  in  an  Atmosphere  of  Humid  Air  and  500 
ppm  Ethylene 


Res.  Time 

(sec) 

NO  Cone. 

(ppm) 

NO2  Cone. 

(ppm) 

Mean  NO 

Mean  NO2 

Runl 

Run2 

Runl 

Run  2 

7.2 

70.3 

74.1 

13.2 

10.4 

72.2 

11.8 

14.5 

68.2 

26.7 

8.5 

15.9 

47.4 

12.2 

29.0 

20.6 

22.2 

7.9 

7.7 

21.4 

7.8 

44.0 

34.6 

13.4 

1.5 

2.4 

24 

1.9 

Table  B.ll:  NO  and  NO2  Data  at  40  kV  in  an  Atmosphere  of  Humid  Air  and  500 
ppm  Ethylene 


Table  B.12:  NO  and  NO2  Data  at  50  kV  in  an  Atmosphere  of  Humid  Air  and  500 
ppm  Ethylene 


Res.  Time 

(sec) 

NO  Cone. 

(ppm) 

NO2  Cone. 

(ppm) 

Mean  NO 

Mean  NO2 

Runl 

Run2 

Runl 

Run  2 

7.2 

15.6 

17.4 

24.2 

24.6 

16.5 

24.4 

14.5 

0 

0 

11.3 

13.8 

0 

12.5 

29.0 

0 

0 

9.8 

14.5 

0 

12.1 

44.0 

0 

0 

2.0 

1.2 

0 

1.6 

APPENDIX  C 


OZONE  FORMATION  DATA 


Table  C.l;  O3  Data  at  30  kV  in  an  Atmosphere  of  Dry  Air 


Res.  Time 

O3  Cone. 

Mean  O3 

(sec) 

(ppm) 

Runl 

Run2 

7.2 

0 

3 

1.5 

14..5 

9 

0 

4.5 

29 

12 

15 

13.5 

44.0 

6 

9 

7.5 

Table  C.2:  O3  Data  at  35  kV  in  an  Atmosphere  of  Dry  Air 


Res.  Time 

O3  Cone. 

Mean  O3 

(sec) 

(PP 

m) 

Runl 

Run2 

7.2 

33 

21 

27 

14.5 

15 

12 

13.5 

29.0 

27 

21 

24 

44.0 

24 

78 

41 

Table  C.3;  O3  Data  at  40  kV  in  an  Atmosphere  of  Dry  Air 


Res.  Time 

O3C 

lone. 

Mean  O3 

(sec) 

(ppm) 

Runl 

Run2 

7.2 

766 

96 

81 

14.5 

45 

81 

63 

29.0 

294 

183 

2.38.5 

44.0 

453 

459 

456 

Table  C.4;  O3  Data  at  50  kV  in  an  Atmosphere  of  Dry  Air 


Res.  Time 

O3  Cone. 

Mean  O3 

(sec) 

(ppm) 

Runl 

Run2 

7.2 

276 

258 

267 

14.5 

375 

363 

369 

29 

1089 

948 

1018.5 

44.0 

1503 

1236 

1369.5 

APPENDIX  D 


COMPUTER  CODES 


D.l  Procedure  of  Computations 

The  pulsed  nature  of  the  discharge  is  simulated  in  the  program  by  considering 
two  reaction  sets  to  calculate  the  concentration  profiles.  During  the  duration  of 
pulsing,  the  reaction  set  which  has  electron-gas  dissociation  reactions  are  used.  The 
time  duration  of  the  pulse  is  assumed  to  be  one  microsecond.  Since  the  frequency  of 
pulsing  is  60  Hz,  the  time  duration  between  the  pulsing  is  16.667  milliseconds.  To 
calculate  the  concentration  profile  for  the  time  in  between  two  pulses,  the  reaction 
set  which  does  not  have  the  electron-gas  reactions  are  considered. 

Initially  the  program  calculates  the  concentration  profile  of  different  species  using 
the  initial  concentrations  given  in  conp.inp  file  and  the  set  of  reactions  which  has 
the  electron-gas  reactions  given  in  chem.inp  for  a  period  of  1  microsecond.  After 
this  time,  since  the  pulsing  is  no  longer  there,  the  reaction  set  is  changed  to  the 
set  which  does  not  have  any  electron-gas  reactions  to  calculate  the  concentration 
profile  for  the  next  16.667  milliseconds.  The  reaction  set  is  again  interchanged  after 
16.667  milliseconds  to  calculate  the  concentration  profile  for  the  next  1  microsec¬ 
ond.  This  procedure  if  repeated  120  times  corresponding  to  double  the  frequency 
of  pulsing  would  give  the  concentration  of  different  species  after  a  gas  residence 
time  of  1  second.  The  concentration  profile  after  each  second  is  copied  into  the  file 
called  cone. out.  The  procedure  is  repeated  till  the  gas  residence  time  in  the  reactor 


is  reached.  The  program  is  repeated  in  this  case  for  7200  times  to  calculate  the 
concentrations  for  a  gas  residence  time  of  60  seconds. 

The  plots  shown  in  the  model  use  the  cone. out  to  show  the  concentrations  of 
different  species  as  a  function  of  residence  time. 

The  program  cycle. sh  is  used  to  repeat  the  iteration  <200  times.  The  program 
conp.sh  is  the  shell  program  which  uses  the  fortran  programs  (conp.f,  vode.f,  cklib.f, 
math.f)  and  the  input  files  to  calculate  the  concentration. 

The  program  conp.f  is  the  actual  fortran  program  which  calls  the  ODE  solver 
vode.f  and  other  subroutines  in  cklib.f  and  math.f  to  solve  the  ODE.  The  ODE 
solver  and  the  other  subroutines  are  a  part  of  the  software  called  CHEMKIN  (Kee, 
et  ah,  1994). 
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D.2  cycle.sh 

#  This  shell  pr-ograns  nuxs  cconp.sh  program  repeatedly  till  the 

#  final  time  correspond*  to  ||he  final  time  in  the  reactor. 

#  In  this  ease,  the  prograii|lt runs  for  60  seconds. 

#  After  that  the  initial  coifeantrat ions  are  added  to  the  cone. out 

#  to  get  the  concentration  j^ofilo, 

#!/bin/csh  -f 

set  icount  -  l 

set  iend  =  7201 

while  ($icount  !=  $iend) 

I? 

sh  conp.sh  conp 


®  acount++ 


D.3  conp.sh 


#  This  shell  program  uses  the  conp.f  program  to  calculate  the  concentration 

#  profile.  It  uses  a  ODE  solver  called  vode.f  and  list  of  subroutines 

#  listed  under  cklib.f  files  from  the  software  called  CHEMKIM. 

#  This  is  shell  program  calculates  the  concentration  profile  till  the 

#  final  time,  interchanges  the  reaction  set  making  the  computed  final 

#  concentration  of  the  initial  concentration  for  the  next  set. 

#  The  outputfile  are  conp.out,  cone. out. 

#  cone. out  has  the  concentration  profile  of  the  species  after  every  second. 

# ! /bin/sh 

#  to  execute:  sh  conp.sh  logname  & 
sh  1>  ${l}.log  2>&1  «  EMDSH 

set  -X 

#cd  /scr/$LOGNAME  #go  to  user's  scratch  directory 

#mkdir  "${!}$$"  #make  subdirectory  /myrun## 

#cd  "${!>$$"  #go  to  /myrun## 

cat  <<  EOF  >  makefile 

include  chemmake.h 

OBJS  =  conp.o  cklib.o  vode.o  math.o 

INPS  =  therm.dat  chem.inp  conp.inp 

OUTS  =  chem.out  chem.bin  conp.out 

EXES  =  chem.exe  conp.exe 

chem.exe:  ckinterp.o 

$(LINK)  chem.exe  ckinterp.o 
conp.exe:  $(0BJS) 

$(LINK)  conp.exe  $(0BJS) 


D.4  conp.f 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


Modeling  of  NO  Removal  From  Various  Gas  Compositions  Using  A 
Pulsed  Streamer  Corona  Discharge  Reactor. 


Input  :  Pressure,  Temperature,  Initial  Mole  Fractions  of 
Different  Species,  Initial  Time,  Final  Time  and 
Interval  in  a  File  Called  conp.inp 
:  Set  of  Chemical  Reactions  Considered  for  the  MOx  Removal 
in  a  File  Called  chem.inp 

Output:  Concentration  Profile  of  Different  Species  till  the 

Final  Time  at  Requested  Intervals  in  a  File  Called  Chem.out 


C  This  Program  calls  the  ODE  solver  to  solve  for  the  given  input 


C  This  program  is  for  the  mole  fractions  as  its  input  . . 

PROGRAM  CONP 
C 

C  Integration  of  adiabatic,  constant  pressure  kinetics  problems 

C 

C  VERSION  1.2: 

C  1 .  Implement  new  VODE  solver 

C 

(;;=it**+*precision  >  double 

IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z),  INTEGER(I-M) 
precision  >  double 
C*****precision  >  single 


c  IMPLICIT  REAL  (A-H,0-Z),  INTEGER  (I-M) 

Q**=it**END  precision  >  single 
C 

PARAMETER  (LENIWK=4000 ,  LENRWK=8000,  LENCWK=500,  MK=5,  NLMAX=55, 

1  LIN=5,  L0UT=6,  LINKCK=25,  KMAX=50,  IT0L=1,  I0PT=0, 

2  RT0L=1.0E-6,  ITASK=1,  AT0L=1 . OE-20) 

C 

DIMENSION  IWORK(LENIWK),  RWORK (LENRWK) ,  X(KMAX) ,  Z(KMAX) 
CHARACTER  CW0RK(LENCWK)*16,  KSYM(KMAX)*16 ,  LINE*80 
LOGICAL  KERR,  lERR 
EXTERNAL  FUN 
C 

COMMON  /RCONS/  P,  RU 
COMMON  /ICONS/  KK,  NWT,  NH,  NWDOT 
C 

DATA  KERR/ . FALSE ./ ,  X/KMAX*0.0/,  KSYM/KMAX* ^  V 

C  Initialize  CHEMKIN 

OPEN ( 19 ,  FILE="Conp . dat " , STATUS="UNKNOWN" ) 

OPEN  (22,FILE="C0NP.DAT1",STATUS="UNKN0WN") 

OPEN  (LINKCK,  FORM=' UNFORMATTED \  STATUS= ' UNKNOWN ' , 

1  FILE='chem.binO 

CALL  CKLEN  (LINKCK,  LOUT,  LENI ,  LENR,  LENC) 

CALL  CKINIT  (LENIWK,  LENRWK,  LENCWK,  LINKCK,  LOUT,  IWORK, 

1  RWORK,  CWORK) 

CLOSE  (LINKCK) 

CALL  CKINDX  (IWORK,  RWORK,  MM,  KK,  II,  NFIT) 


NEQ  =  KK  +  1 

LRW  =  22  +  9*NEQ  +  2*NEQ**2 
MVODE  =  LENR  +  1 
NWT  =  NVODE  +  LRW 
NH  =  NWT  +  KK 
NWDOT  =  NH  +  KK 
NTOT  =  NWDOT+  KK  -  1 
LIW  =  30  +  NEQ 
IVODE  =  LENT  +  1 
ITOT  =  IVODE  +  LIW  -  1 
IF  (KK  .GT.  KMAX)  THEN 
WRITE  (LOUT,  *) 

1  '  Error. . .KMAX  too  small. . .must  be  at  least  ’ ,KK 

KERR  =  .TRUE. 

ENDIF 

IF  (LEMRWK  .LT.  NTOT)  THEN 
KERR  =  .TRUE. 

WRITE  (LOUT,  *) 

1  '  Error. . .LENRWK  too  small... must  be  at  least',  NTOT 

ENDIF 

IF  (LENIWK  .LT.  ITOT)  THEN 
KERR  =  .TRUE. 

WRITE  (LOUT,  *) 

1  '  Error. . .LENIWK  too  small... must  be  at  least',  ITOT 

ENDIF 

IF  (KERR)  STOP 

CALL  CKSYMS  (CWORK,  LOUT,  KSYM,  lERR) 


IF  (lERR)  KERR  =  .TRUE. 

CALL  CKWT  (IWORK,  RWORK,  RWORK(NWT)) 

CALL  CKRP  (IWORK,  RWORK,  RU,  RUC,  PATH) 

READ  (LIN,  ♦)  PA,  T 
WRITE (19,*)  PA,  T 
P  =  PA*PATM 

C  Initial  non-zero  moles 

40  CONTINUE 
LINE  =  ^  ' 

READ  (LIN,  ' (A) '  ,  END=45)  LINE 
ILEN  =  INDEX  (LINE,  ' ! ') 

IF  (ILEN  .EQ.  1)  GO  TO  40 
ILEN  =  ILEN  -  1 

IF  (ILEN  .LE.  0)  ILEN  =  LEN(LINE) 

IF  (IMDEX(LINE(:ILEN) ,  'END')  .EQ.  0)  THEN 
IF  (LINE(;ILEN)  .NE.  '  ')  THEN 

CALL  CKSNUM  (LINE (: ILEN) ,  1,  LOUT,  KSYM,  KK,  KNUM, 
1  NVAL,  VAL,  lERR) 

IF  (lERR)  THEN 

WRITE  (LOUT,*)  '  Error  reading  moles...' 

KERR  =  .TRUE. 

ELSE 

X(KNUM)  =  VAL 
ENDIF 
ENDIF 
GO  TO  40 


ENDIF 


45  CONTINUE 


C  Final  time  and  print  interval 

READ  (LIN,  *)  TT1,T2,  DT,  III 
WRITE  (22,7100)  (KSYM(K),  K=1,KK) 

IF  (KERR)  STOP 

C  Initial  conditions  and  mass  fractions 

Z(l)  =  T 

CALL  CKXTY  (X,  IWORK,  RWORK,  Z(2)) 

C  Integration  control  parameters  for  VODE 

TT2  =  TTl 
MF  =  22 
ISTATE=  1 
NLINES=NLMAX  +  1 
C  Integration  loop 

250  CONTINUE 
T  =  Z(l) 

CALL  CKYTX  (Z(2),  IWORK,  RWORK,  X) 

CALL  CKMMWY(Z(2),  IWORK,  RWORK,  WW) 

WRITE  (LOUT,  7105)  TTl,  T,  (X(K),  K=1,KK),  WW 

IF  (TT2  .GE.  T2)  THEN 

WRITE (19, 70 12) (KSYM(K) ,X(K) ,  K=1,KK) 

WRITE (19,*) 'END' 

JJJ=III/2 

QQQ=JJJ*2 

IF(Qqq.EQ.III)  THEN 
XINT=lE-6 


T4=TT2 


TT1=T2 


T2=T2+XINT 

T3=XINT/10 

IF  (M0D(III,120) .EQ.O)  THEN 
WRITE(25,7105)  TT2,T,(X(K),  K=  1,KK),  WW 
END  IF 
ELSE 

XINT=16.667E-3 

TT1=T2 

T2=T2+XINT 

T3=XINT/10 

END  IF 

III=III+1 

WRITE(19,*)  TTl,  T2,T3,  III 

STOP 

ELSE 

TT2  =  MIN(TT2  +  DT,  T2) 

END  IF 

C  Call  the  differential  equation  solver 
350  CONTINUE 
C*****precision  >  single 
C  CALL  SVODE 

C*****END  precision  >  single 
C*****precision  >  double 
CALL  DVODE 

C*****END  precision  >  double 
* 


(FUN,  NEQ,  Z,  TTl,  TT2,  ITOL,  RTOL,  ATOL,  ITASK, 


1  ISTATE,  lOPT,  RWORK(NVODE),  LRW,  IWORK(IVODE) , 

2  LIW,  JAC,  MF,  RWORK,  IWORK) 

IF  (ISTATE  .LE.  -2)  THEN 

IF  (ISTATE  .EQ.  -1)  THEN 
ISTATE  =  2 
GO  TO  350 
ELSE 

WRITE  (LOUT,*)  '  ISTATE=' , ISTATE 
STOP 
ENDIF 
END  IF 
GO  TO  250 
C  FORMATS 

7003  FORMAT  (IHl) 

7100  FORMAT  (3X,  'T(SEC) ’ ,  9X,  'TMP(K)',  9X,  20(4X,A10)) 

7105  FORMAT  (50E16.8) 

7110  FORMAT  (26X,  5(1X,A10)) 

7115  FORMAT  (22X,  10E11.3) 

7012  FORMAT  ((1X,A10) ,4X,1E14.8) 

END 

SUBROUTINE  FUN  (N,  TIME,  Z,  ZP,  RPAR,  IPAR) 

C*****precision  >  double 

IMPLICIT  DOUBLE  PRECISION (A-H, 0-Z) ,  INTEGER(I-N) 

C*****END  precision  >  double 
C*****precision  >  single 

C  IMPLICIT  REAL  (A-H, 0-Z),  INTEGER(I-N) 

C*****END  precision  >  single 


COMMON  /RCONS/  P,  RU 

COMMON  /ICONS/  KK,  NWT,  NH,  NWDOT 

DIMENSION  Z(*)  ,  ZP(*),  RPAR(*) ,  IPAR(*) 

C 

C  Variables  in  Z  are:  Z(l)  =  T 

C  Z(K+1)  =  Y(K) 

C 

c  Call  CHEMKIN  subroutines 

C 

CALL  CKRHOY  (P,  Z(l),  Z(2),  IPAR,  RPAR,  RHO) 

CALL  CKCPBS  (Z(l)  ,  Z(2) ,  IPAR,  RPAR,  CPB) 

CALL  CKWYP  (P,  Z(l),  Z(2),  IPAR,  RPAR,  RPAR(NWDOT)) 
CALL  CKHMS  (Z(l),  IPAR,  RPAR,  RPAR(NH)) 

C 

C  Form  governing  equation 

C 

SUM  =  0.0 
DO  100  K  =  1,  KK 

H  =  RPAR(NH  +  K  -  1) 

WDOT  =  RPAR (NWDOT  +  K  -  1) 

WT  =  RPAR(NWT  +  K  -  1) 

ZP(K+1)  =  WDOT  *  WT  /  RHO 

C  SUM  =  SUM  +  H  *  WDOT  ♦  WT 

100  CONTINUE 
RETURN 


END 


D.5  conp.inp 


1.00000000000000 


02 
0 
N 

N2 
MO 
03 
N02 
NOS 
N20 
END 

0  lE-06  lE-07  1 


0.2099 

0. 

0. 

0.79 

0.0001 

0. 

0. 

0. 

0. 


298.000000000000 


